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We present exact calculations of the partition function Z of the g-state 
Potts model and its generalization to real q, the random cluster model, for 
arbitrary temperature on n-vertex ladder graphs with free, cyclic, and Mobius 
longitudinal boundary conditions. These partition functions are equivalent to 
Tutte/Whitney polynomials for these graphs. The free energy is calculated 
exactly for the infinite-length limit of these ladder graphs and the thermo- 
dynamics is discussed. By comparison with strip graphs of other widths, we 
analyze how the singularities at the zero-temperature critical point of the 
ferromagnet on infinite-length, finite-width strips depend on the width. We 
point out and study the following noncommutativity at certain special val- 
ues Qs'- lim„^oo hmg^g^ Z^/" ^ lim^^g^ lim„^oo It is shown that the 
Potts/random cluster antiferromagnet on both the infinite-length line and 
ladder graphs with cyclic or Mobius boundary conditions exhibits a phase 
transition at finite temperature if < g < 2, but with unphysical properties, 
including negative specific heat and non-existence, in the low-temperature 
phase, of an n — > oo limit for thermodynamic functions that is independent of 
boundary conditions. Considering the full generalization to arbitrary complex 
q and temperature, we determine the singular locus B in the corresponding 
space, arising as the accumulation set of partition function zeros as n — > co. In 
particular, we study the connection with the T = limit of the Potts antifer- 
romagnet where B reduces to the accumulation set of chromatic zeros. Certain 
properties of the complex-temperature phase diagrams are shown to exhibit 
close connections with those of the model on the square lattice, showing that 
exact solutions on infinite-length strips provide a way of gaining insight into 
these complex-temperature phase diagrams. 
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I. INTRODUCTION 



The g-state Potts model has served as a valuable model for the study of phase transitions 
and critical phenomena [|l],@. On a lattice, or, more generally, on a graph G, at temperature 
T, this model is defined by the partition function 

Z(G,g,t;) = 5^e-^^ (1.1) 
with the (zero-field) Hamiltonian 

n = -jY,5„.., (1.2) 

where (Tj = 1, g are the spin variables on each vertex i E G; (3 = {kBT)~^; and (zj) denotes 
pairs of adjacent vertices. The graph G = G(y, E) is defined by its vertex set V and its edge 
set E] we denote the number of vertices of G as = n{G) = \V\ and the number of edges of 
G as e(G) = We use the notation 

K = (3J (1.3) 



a = u-' = (1.4) 

and 

v = a-l (1.5) 

so that the physical ranges are (i) a > 1, i.e., f > corresponding to oo > T > for the 
Potts ferromagnet, and (ii) < a < 1, i.e., —1 < f < 0, corresponding to < T < oo for 
the Potts antiferromagnet. An equivalent expression for Z is 

ZiG,q,v) = J2U(^ + vS..,^,) ■ (1-6) 

One defines the (reduced) free energy per site / = — where F is the actual free energy, 
via 

f{{G},q,v)= hm ln[Z(G,g,t;)V"] . (1.7) 

n^oo 

Let G' = {V, E') be a spanning subgraph of G, i.e. a subgraph having the same vertex 
set V and an edge set E' C E. Then Z{G, q, v) can be written as the sum 0- |^ 
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Z{G,q,v)= ^gMG')^e(G') (1^^ 



G'CG 

where k{G') denotes the number of connected components of G' . The formula ( p. .81) enables 
one to generalize q from Z+ to IR+ (keeping f in its physical range). This generalization 
is the random cluster model 0. The formula ( |1.8|) shows that Z{G,q,v) is a polyno- 
mial in q and f (equivalently, a) with maximum degrees max{(ie5'g(Z(G, g, f ))} = n{G) 
and max{cie(y'^(Z(G', g, t>))} = e{G). The minumum degrees are mm{degq{Z{G,q,v))} = 
k{G), which is equal to 1 for the graphs of interest here (since they are connected), and 
mm{degjj{Z{G,q,v))} = 0, so 

n{G) e(G) 

Z{G,q,v)= J2 E^-^'^' (1-9) 

r=k{G) s=0 

with Zrs > 0. 

The Potts model partition function on a graph G is essentially equivalent to the Tutte 
polynomial j^- |TT] and Whitney rank polynomial |Q, 0, ||T^- \T^\ for this graph, as dis- 



cussed in the appendix. As a consequence, there are many interesting connections between 
properties of this partition function and various graph-theoretic quantities. 

The Potts model has never been solved exactly for arbitrary temperature on lattices of 
dimensionality d > 2 except for special d = 2, q = 2 case in which it is equivalent to the 
solvable 2D Ising model []15| (with the redefinition Jpotts = 2J/smg)- Knowledge about the 



Potts model includes exact calculations of the critical exponents and critical value of the free 
energy for the 2D ferromagnet for the range g < 4 where it has a second-order transition; 
conformal algebra properties for the same range of q; the latent heat at the transition point 
for g > 5; and certain formulas for the critical point P,p^. There is thus motivation for 
studies that can give further insight into the Potts model. Among these are exact results for 
the partition function and free energy that one can obtain for infinite-length, finite-width 
strips with various boundary conditions. We shall present such results in this paper. 



One of the interesting features of the Potts model is that the antiferromagnet (AF) 
exhibits nonzero ground state entropy (without frustration) for sufficiently large g on a 
given lattice or graph G, and serves as a valuable model for the study of this phenomenon. 
The phenomenon of nonzero ground state entropy, 5*0 > 0, is an exception to the third law 
of thermodynamics [|1^|18|. This is equivalent to a ground state degeneracy per site (vertex), 
W > 1, since 5*0 = ks In W. The zero-temperature partition function of the above-mentioned 
g-state Potts antiferromagnet (PAF) on G satisfies 
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Z(G, g, T = 0)pAF = Z{G, q,v = -1) = P{G, q) 



(1.10) 



where P(G, q) is the chromatic polynomial (in q) expressing the number of ways of coloring 
the vertices of the graph G with q colors such that no two adjacent vertices have the same 



color ||3| JT9|j20[] . The minimum (integral) number of colors necessary for this coloring is the 



chromatic number of G, denoted x{G)- Thus 

iy({G},g) = lim P(G,g)^/" (1.11) 

where we use the symbol {G} to denote lim^^oo G for a given family of graphs. 

Since Z{G^ q, v) is a polynomial in q and f , or equivalently, a, one can generalize q from 
Z_|_ not just to M-i. but to C and a from its physical ferromagnetic and antiferromagnetic 
ranges 1 < a < oo and 0<a<ltoaGC A subset of the zeros of Z in the two-complex 
dimensional space defined by the pair of variables (g, a) can form an accumulation set in 
the n — » oo limit, denoted which is the continuous locus of points where the free energy 
is nonanalytic. As will be discussed below, this locus is determined as the solution to a 
certain {G}-dependent equation. For a given value of a, one can consider this locus in the 
q plane, and we denote it as Bq{{G},a). In the special case a = (f = —1) where the 
partition function is equal to the chromatic polynomial, the zeros in q are the chromatic 
zeros, and i3q({G},a = 0) is their continuous accumulation set in the n — oo limit P^ - 



P7| ; we have determined these accumulation sets exactly for various families of graphs in 
a series of papers. Other properties of chromatic zeros such as zero-free regions for general 
graphs, are of mathematical interest (see, e.g., |^,^,|^), although we shall not focus on them 
here. For a given value of g, we shall study the continuous accumulation set of the zeros 
of Z{G,q,v) in the a plane; this will be denoted Ba{{G},q). It will often be convenient to 
consider the equivalent locus in the u = 1/a plane, namely Ba{{G},q). We shall sometimes 
write Bq{{G}, a) simply as Bg when {G} and q are clear from the context, and similarly with 
Ba and Ba- A subtlety in the definition of this locus will be discussed in the next section. 

One gains a unified understanding of the separate loci Bq{{G}) for fixed a and Ba{{G}) 
for fixed q by relating these as different slices of the locus B in the space defined by (g, a). 
This is similar to the insight that one gained in studies of accumulation sets of zeros of the 
partition functions of the Ising model in the space defined by (a, /x), where fi = e'"^^^ 
with H being the external field p5|J5m , which generalized the Yang-Lee zeros (zeros in for 



fixed physical a) |^ and Fisher zeros (zeros in a for fixed often H = 0) 

In our earlier works on Bg{{G}) for a = 0, we had denoted the maximal region in the 
complex g plane to which one can analytically continue the function iy({G}, g) from physical 
values where there is nonzero ground state entropy as -Ri. The maximal value of g where 
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B intersects the (positive) real axis was labelled qc{{G}). Thus, region Ri includes the 
positive real axis for q > qc{{G}). Correspondingly, in our works on complex-temperature 
properties of spin models, we had labelled the complex-temperature extension (CTE) of 
the physical paramagnetic phase as (CTE)PM, which will simply be denoted PM here, the 
extension being understood, and similarly with ferromagnetic (FM) and antiferromagnetic 
( AFM) ; other complex-temperature phases, having no overlap with any physical phase, were 
denoted Oj (for "other"), with j indexing the particular phase [Q. Here we shall continue 
to use this notation for the respective slices of B in the q and a or m planes. 



In this paper we shall present exact calculations of the Potts/random cluster partition 
function for strips of the square lattice with arbitrary length L^. and width Ly = 2, i.e ladder 
graphs, having free, periodic (= cyclic), and Mobius longitudinal (x-direction) boundary 
conditions. These famihes of graphs are denoted, respectively, as Sm (for open strip), 
(for ladder) , and ML^ (for Mobius ladder) , where L^. = m + 1 for Sm (following our labelling 
convention in |^) and = m for Lm and MLm- One has n{Sm) = 2(m + 2) and n{Lm) = 



n{MLm) = 2m. It will also be instructive to use the well-known solutions for the partition 
function on the tree and circuit graphs to illustrate some points. We shall discuss several 
items and investigate several questions about the Potts/random cluster model: 

1. We analyze the thermodynamic behavior of the Potts model (for g > 2) on the infinite- 
length, Ly = 2 strip and compare it with the known behavior on the line. In particular, 
we study the zero-temperature critical point of the Potts ferromagnet and discuss how 
the critical singularities (which are essential singularities in temperature) depend on 
q and the width of the strip. For reference, one may recall that as one part of his 
original paper, Onsager used his closed-form solution to the partition function of the 



Ising model to study it for Ly x oo strips of the square lattice [jl5l. The difference, of 
course, is that for the Potts model with q ^ 2, one does not have a general closed-form 
solution for Z on a Ly x L^ grid with arbitrary Ly and L^; indeed if one did, one would 
have solved the model on the square lattice. 

2. We shall show that the Potts/random cluster model with < g < 2 on the n oo 
limit of the circuit, ladder, and Mobius ladder graphs exhibits a finite-temperature 
phase transition but with unphysical properties in the low-temperature phase, including 
negative specific heat, negative partition function, and non-existence of an n ^ oo limit 
that is independent of boundary conditions. 

3. We shall point out and illustrate a certain noncommutativity in the definition of the 
free energy for the Potts/random cluster model. 
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4. For a given family of graphs G and its n — >^ oo limit, {G}, we shall explore the 
nature of the nonanalyticities of the free energy in q and the temperature variable u. 
Some questions related to this are: what is the locus Bq for various values of u and 
the locus Bu for various values of q (and the loci {Bu)qn and {Bu)nq)i for the special 
values q = qs where these differ)? The strip graphs that we consider here are useful 
for this study since they are wide enough to exhibit a number of important features 
but narrow enough so that the terms, denoted Xj, whose powers contribute to the 



partition function for a given family of graphs (see eq. ( 2.18 )) can be calculated as 



explicit algebraic functions. As our previous studies of asymptotic limits of chromatic 



polynomials on various infinite-length, finite- width strips have shown ||3^j3^j41|j4t^ , as 
the width of the strip increases, one encounters algebraic equations defining the Aj's 
that increase in degree so that these Xj can involve cube roots, fourth roots, and, for 
equations higher than fourth degree, one cannot obtain exact analytic expressions for 
them, making it more cumbersome to calculate the locus B, which is the solution to 
the degeneracy in magnitude of different dominant Aj's. 

5. Starting from our previous determination of Bq{{G}) for the zero-temperature limit 
of the Potts antiferromagnet, we explore how this locus changes as on increases T in 
the range < T < oo. In cases where this locus separates the q plane into different 
regions for T = 0, does it continue to do so? How does the point qc{{G}) vary with 
temperature? Since we are now dealing with a singular locus in C^, we can investigate 
how the various features of the slice of B{{G}) in the q plane relate to features of the 
slice B{{G}) in the plane of the temperature variable, u (or the equivalent variable a). 

6. We shall show that certain features of the complex-temperature phase diagrams of the 
infinite-length, finite-width strip graphs considered here can give insight into analogous 
features of the Potts model on the square lattice. 

7. Just as was true for our earlier studies of chromatic polynomials and their asymptotic 
limits, it is of interest to explore the effects of different boundary conditions on the 
singular locus B{{G}), and we do this. 

8. For mathematically inclined readers, we shall give the Tutte polynomials that are 
equivalent to the Potts model partition functions that we have calculated for the cyclic 
and Mobius strip graphs and extract special values that are of specific graph-theoretic 
interest. 



Conference reports on some of the material in this paper were given in ||43|J4^ . In col- 
laboration with H. Kluepfel, we have also carried out calculations of Potts/random cluster 
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partition functions for general T and q on finite patches of several 2D lattices (see also 
53|); our work here complements these calculations on finite patches in that we obtain ex- 



act results for strip graphs of arbitrarily great length, and the nonanalyticities in the limit 
n — s> oo. 



II. SOME GENERAL CONSIDERATIONS 
A. Basic Properties of Z 

For our later analysis, it will be useful to record some basic properties of the Potts model 
partition function. Assuming g > 0, one observes that for the Potts ferromagnet, since 
f > 0, each term in the sum ( |1.8| ) is positive, and consequently, Z{G, q, v) does not have any 
zeros on the positive real q axis for the physical temperature range. On higher-dimensional 
lattices where the ferromagnet has a finite-temperature phase transition, zeros will coalesce 
and pinch the real positive q axis to form a region boundary, but this does not happen in 
the ID case and the infinite-length, finite-width strips, which are quasi-lD systems. 

One may ask what factors Z has in general. From eq. ( |1.8| ) it follows that for an arbitrary 
graph G, 

Z{G,q = 0,v) = (2.1) 

and since Z{G, q, v) is a polynomial, this implies that Z{G, q, v) always has an overall factor 
of q. For the families of graphs studied here, this is, in general, the only overall factor that 
Z(G,q,v) has. For the special case v = —1, the resultant chromatic polynomial Z{G,q,v = 
— 1) = P{G, q) has the additional factors nsi'i''~^(5' ~ Another general result is that 

Z(G, q = l,v)= v<^''^ = a"(^) . (2.2) 

G'CG 

For temperature T = oo, i.e., f = 0, we have 

Z{G,q,v = Q) = q<^^ . (2.3) 

For the Ising case g = 2, if G is bipartite, 

Z{Gup., q = 2,a) = a^^"""^ Z^Gup., q = 2, 1/a) (2.4) 

where we have written Z here as a function of a. This is the well-known equivalence of the 
Ising ferromagnet (F) and antiferromagnet (AF) on a bipartite lattice, when one makes the 
replacement J — J, and hence K —K, a ^ 1/a. As a consequence, the zeros of Z in a 
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for q = 2 are invariant under the inversion mapping a ^ 1/a. Among the famihes of graphs 
considered here, the following are bipartite (equivalently, have chromatic number x = 2): 
tree graph Tm for any m; circuit graph Cm and cyclic ladder graph L„i with Ly = 2 for even 
m; and Mobius ladder graph MLm for odd m. In contrast, the cyclic ladder graph Lm for 
odd m and the Mobius ladder graph MLm for even m have x = 3. 

For a strip of the square lattice with width Ly and cyclic or Mobius longitudinal boundary 
conditions, the average coordination number (degree in the graph-theoretic terminology) 
21inin(G)^oo e(G')/n(G') is 

Aave. = 4 - ^ . (2.5) 

Ly 

For the corresponding strip of width Ly and free boundary conditions, this formula for A 
also holds in the L^ ^ oo limit. Another consequence of the symmetry ( |2.4| ) is that for the 
Ising case the internal energy U and specific heat C satisfy the relations 

[/(Gmp, q = 2, J)f = U{Gmp., q,J^ -J)af - (2.6) 

and 

C{Ghip, q = 2, J)p = C{Gbip,, q, J ^ —J)af (2.7) 

(where we have taken the n ^ oo limit, in which these results are independent of the 
boundary conditions for physical values of temperature). 

Another basic property, evident from eq. (|1.8| ), is that (i) the zeros of Z{G, q, v) in q for 
real v and hence also the continuous accumulation set Bq are invariant under the complex 
conjugation q q*; (ii) the zeros of Z{G, q, v) in v or equivalently a for real q and hence also 
the continuous accumulation set Ba are invariant under the complex conjugation a — >■ a*. 



B. Noncommutativity in the Random Cluster Model 



Just as we showed the importance of noncommutative limits in our earlier work on chro- 
matic polynomials (eq. (1.9) in Ref. |^), so also we encounter an analogous noncommuta- 
tivity here for the general partition function ( |1.8|) of the random cluster model: at certain 
special points qs (typically qs = 0, 1..., x{G)) one has 

lim lim Z{G,q,vy/'^ ^ lim lim Z(G,g,t;)^/" . (2.8) 

n— >oo q^qs q^qs n^oo 

Clearly, no such issue of noncommutativity arises if one restricts to positive integer q values 
and uses the Potts model definition (|1.1| ), ( p.. 21) . However, for the general random cluster 
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model, whenever Z{G,q,v) has a factor (q — qs)^" with finite multiphcity fig (where typically 
/is = 1 here), one encounters this noncommutativity. It can also occur when such a factor 
appears as a coefficient of one of the terms (Agj)™' contributing to Z{G, q, v) (see eq. ( p.l8|) 



below). The reason for this noncommutativity is analogous to that which we discussed earlier 
in the special case (a = 0) of the chromatic polynomial ||2l]]; it is a consequence of the basic 
result 

hm (g - g,)'^^/" = | J ^ . (2.9) 

We shall illustrate this with our exact results to be discussed below. Because of this non- 
commutativity, the formal definition ( |1.7|) is, in general, insufficient to define the free energy 
/ at these special points g^; it is necessary to specify the order of the limits that one uses in 
eq. ( |2.8| ). We denote the two definitions using different orders of limits as fgn and fnq- 

fnci{G},q,v)= lim hm n-'\nZiG,q,v) (2.10) 

n— >oo q^qs 



fqn{{G}, g, v) = lim lim n'^ In Z{G, g, v) . (2.11) 

For the zero-temperature Potts/random cluster antiferromagnet case a = (f = —1), the 
same ordering ambiguity affects the formal equation (|1 . 1 1| ) . In Ref. I^Tj and our subsequent 



works on chromatic polynomials and the above-mentioned zero-temperature antiferromag- 
netic limit, it was convenient to use the ordering W{{G},qs) = limg^^^ lim^^oo P(G, g)^/" 
since this avoids certain discontinuities in W that would be present with the opposite or- 
der of limits. In the present work on the full temperature-dependent random cluster model 
partition function, we shall consider both orders of limits and comment on the differences 
where appropriate. Of course in discussions of the usual g-state Potts model (with positive 
integer g), one automatically uses the definition in eq. (|1.1|) with ( p.. 21) and no issue of orders 
of limits arises, as it does in the random cluster model with real g. 

As a consequence of the noncommutativity ( p.8| ), it follows that for the special set of 
points q = qs one must distinguish between (i) {Ba{{G}, qs))nq, the continuous accumulation 
set of the zeros of Z{G, g, v) obtained by first setting q = qs and then taking n oo, and 
(n) {Ba{{G},qs))qn, the continuous accumulation set of the zeros of Z{G,q,v) obtained by 
first taking n —>■ oo, and then taking g — > g^. For these special points, 

{Ba{{G}, qs))nq + (Sa({G}, g.)),. . (2.12) 

From eq. (p.l|), it follows that for any G, 
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exp(/„,) = for g = (2.13) 

and thus 

{BaU = ^ for g = 0. (2.14) 

However, for many families of graphs, including the circuit graph C„, and cyclic and Mobius 
strips Lm and MLm, if we take n — oo first and then g — 0, we find that {Ba)qn is nontrivial. 
For these families of graphs, with this order of limits, although the free energy is nonanalytic 
at g = 0, it is continuous, and we find that, in general, exp(/q„) 7^ at g = 0. Similarly, 
from (|2.2|) we have, for any G, 

{Ba)nq = ^ for g = l (2.15) 

since all of the zeros of Z occur at the single discrete point a = (and in the case of a graph 
G with no edges, Z = 1 with no zeros). However, as the simple case of the circuit graph 
below will show, {Ba)qn is, in general, nontrivial. 

We shall also find that {Ba)nq 7^ {Ba)qn for g = 2 for the infinite-length, Ly = 2 width 
strip graphs {L} and {ML}. As stated, this noncommutativity can, in general, occur at 
integer values of g up to and including g = x{G). However, although x = 3 for and 
Lm with odd m and for ML^, with even m, there is no noncommutativity at g = 3 in these 
cases. This can be seen as a consequence of the fact that one can take the limit m — >^ 00 on 
even values of m. 

In the g = 2 Ising case, as a consequence of the relation ( p.4|) , the locus {Ba)nq is invariant 
under the inversion map a —>■ 1/a for the n 00 limit of a sequence of bipartite graphs of 
type G: 

iBa)nq{{Gbip.}) is invariant under a — > - if g = 2 (2-16) 

ct 

where {Gfejp.} means that for the family of graphs of type G, one can take the limit n ^ 00 
with a sequence of bipartite members of the family G. (For example, for the circuit graphs 
Cn, one can do this by taking n — 00 on even values, and so forth for other families.) As our 
explicit results for the strips {L} and {ML} below will show, the locus obtained with the 
opposite order of limits, {Ba)nq{{Gbip.}), does not, in general, have this inversion symmetry, 
even if g = 2. 

Concerning the cases where the continuous locus B may be the nullset 0, some examples 
from chromatic polynomials may be useful. For the complete graph on n vertices (defined 
as the graph in which each point is connected to every other point with edges), if a = 0, 
then Z{Kn, q,v = —1) = P{Kn, q) = 11^=0 (•? ~ ■^)' ^^^^ as n — 00, there is no continous 
accumulation set of the chromatic zeros, so ;Bg = . Another example is provided by the 
strip Sm (see eq. (|5.27| ) below). 
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C. Definition of Free Energy for Complex q and K 



Another matter concerns the definition of / away from physical values of q and where 
Z{G, q, a) can be negative or complex. In these ranges of q and K, there is no canonical 
choice of which 1/n'th root, i.e., which value of r, to pick in eq. ( p. 10 ) or ( p. 11 ): 



Z(G,g,a)^/" = {|Z(G,g,a)|^/V'^+2-*")/"} , r = 0, 1, n - 1 . (2.17) 

where = arg(Z). Thus, we start by considering the free energy / defined for sufficiently 
large (physical) T and integer q > 2 and define the maximal region in the q plane for 
fixed a or in the a plane for fixed q that can be reached by analytic continuation of this 
function. As noted, this is labelled the region Ri in the q plane and the PM phase (and its 
complex-temperature extension) in the u plane. In these regions, the canonical phase choice 
in (|2.17|) is clearly that given by r = 0. This would also be true in physical low-temperature 
broken-symmetry phases such as ferromagnetic (FM) or antiferromagnetic (AFM) phases 
and their complex-temperature extensions, as discussed in [0]; however, such phases do not 
occur in the ID and quasi-lD strip graphs considered here. However, in general, in complex- 
temperature O phases in the u plane and Rj regions with j 7^ 1 in the q plane, only the 
quantity |e-^| = lim^^oo \Z{G, q, a)\^/^ can be determined unambiguously. 



D. General Form of Z for Recursively Defined Graphs 



We find that a general form for the Potts model partition function for the strip graphs 
considered here, or more generally, for recursively defined families of graphs comprised of m 
repeated subunits (e.g. the columns of squares of height Ly vertices that are repeated 
times to form an Lx x Ly strip of a regular lattice with some specified boundary conditions), 
is 

Z{G,q,v) = J2cGA^GAl^^)r (2-18) 

where N\ depends on G. The formula ( p.l8| ) can be understood from the fact that for the 
cyclic case, for q G Z+, Z{G, q, v) can be expressed as the trace of a transfer matrix T: 

Z{G,q,v)=Tr{rn . (2.19) 

Having obtained Z{G, q, v) in this manner, one can then generalize q from Z_|_ to M^. For a 
strip of a regular lattice, this transfer matrix has dimensions M x A/", where M denotes the 
number of possible spin configurations along a transverse slice of the strip; for example, for 
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the square strips of interest here (or for triangular strips, written in the form of a square 
strip with additional diagonal edges), M = q^y . Clearly, T is a symmetric matrix each of 
whose elements is either 1 or a (positive) power of a = f + 1. For physical temperature, 
for which a > 0, T is real and hence can be diagonalized by an orthogonal transformation 
O, yielding Af eigenvalues (some of which may coincide). Generically, the multiplicity of 
a given eigenvalue Xgj, which yields the coefficient Cj for these cyclic graphs, depends on 
q but is independent of v. The result ( p.l8|) applies for both free and periodic transverse 
boundary conditions, given that one uses periodic (cychc) longitudinal boundary conditions. 
Similar arguments based on the transfer matrix yield this result if one uses free transverse 
and Mobius longitudinal, and periodic transverse and Mobius longitudinal (i.e. Klein bottle) 
boundary conditions. 

For strips with open longitudinal boundary conditions, Z{G,q,v) is not a trace, but 
instead, is given by 

Z{G,q,v) = J2{^^\'^''^\^f) (2-20) 

where 5"j and df denote the states of the Ly spins on the initial (i) and final (f) transverse 
edges of the strip. As our explicit solution (given below) for the open ladder graph Sm shows, 
here Cj can depend on both q and v. 

For recursively defined families of graphs G, the result ( |2.18| ) is a generalization to the 
case of the Potts/random cluster model partition function Z{G, q, v) of the Beraha-Kahane- 
Weiss result that the chromatic polynomial P{G,q) = Z{G,q,v = —1) can be written 
in the form 

P{G,q) = J2cp,GA^P,G,r- (2.21) 

Since P{G, q) is a special case of Z{G, q), it follows that 

Nx,p < Nx . (2.22) 

For example, in two well-known cases, (i) tree graph: A^a = A^a,p = 1, (ii) circuit graph: 
Nx = Nx,p = 2. We find (iii) for the open Ly = 2 ladder Sm, Nx = 2, while Nx,p = 



1; (iv) f orthecydic and Mobius strips, Nx = Q while p2[ Nx^p = 4. For the subset of the 
Agj's that remain for P{G,q), the coefficient functions cp^cj = CG,j{v = —1); as remarked 
above, for the cyclic case and our Mobius strip, cqj are independent of v. 

As m — ^> oo, for a given point (g, v) in the space of variables, one Xj will dominate this 
sum; we denote this as the "leading term" A^, where i stands for leading. As one moves to 
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another point (g', f '), it may happen that there is a change in the dominant A, from A^ to, say, 
A^. Consequently, there is a nonanalytic change in the free energy / as it switches from being 
determined by the first dominant A^ to being determined by A^. Thus, the equation for the 
continuous nonanalytic locus B across which / is nonanalytic, is (with the {G} dependence 
indicated explicitly) 

B{{G}) : |Ag,.| = \X'gA ■ (2.23) 

Although / is nonanalytic across B, a consequence of eq. (|2.23|) is that |e-^| is continuous 
across this locus. This is the generalization of the analogous phenomenon for the asymptotic 
limit of chromatic polynomials, or equivalently, the T = limit of the Potts antiferromagnet 
p5| , |26|p8| , |2T|| , i.e. the property that W{{G},q) is nonanalytic (although its magnitude is 
continuous) across Bq{{G}). As noted above, the locus B forms as the continuous accumu- 
lation set of the zeros of Z{G, q, v) in the 2-complex dimensional space (g, v) as n{G) oo. 
This again generalizes the earlier analysis for chromatic polynomials, where Bq formed as the 
continuous accumulation set of the zeros of the chromatic polynomial in the single complex 
variable q. 

It is straightforward to generalize the transfer matrix formalism and hence eqs. (|2.19|) and 



( p.20| ) to the case of the Potts model in an external magnetic field H, where the Hamiltonian 
is 7"^ = —JJ2(ij) '^o-io-j ~ H J2i '^0-1,0-0 (with (To chosen, say, as 1). Hence, the full generalization 
for recursively defined families of graphs is, with r] = e^^ , 

ZiG,q,v,v) = J2^o,MGAl,^,v)r ■ (2.24) 
i=i 

One could proceed to study the singular locus B in the space of the variables {q,v,ri) 
which forms as the continuous accumulation set of the zeros of Z, and the various slices in 
the q, a, and t] planes. Here we restrict ourselves to the zero- field case, rj = 1. 

Typically, for a given point a G C, there will be an infinite set of points in the q plane lying 
on B, and for a given point g € C, there will be an infinite set of points in the a plane lying 
on B. Again, usually (we have commented on some exceptions above, where accumulation 
sets of zeros are discrete), we find that Bg and Ba form curves (and possible line segments) 
in the respective q and x = a or -u planes. This follows from the property that B is the 
solution of an algebraic equation expressing the degeneracy of the leading A's contributing 
to Z. (For higher dimensional lattices the equations defining B can be transcendental instead 
of algebraic.) The loci Bq and Bx, x = a or u, may be connected, or may consist of several 
disconnected components; our results illustrate both types of behavior. 
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The Potts ferromagnet has a zero-temperature phase transition in the ^ oo hmit of 
the strip graphs considered here, and this has the consequence that for cychc and Mobius 
boundary conditions, B passes through the T = point u = 0. It follows that B is non- 
compact in the a plane. Hence, it is usually more convenient to study the slice of B in the 
u = 1/a plane rather than the a plane. In the Ising case q = 2, Ba = Bu and so both are 
noncompact. For the ferromagnet, since a — oo as T — and Z diverges like a'^*-'^^ in this 
limit, we shall use the reduced partition function Zr defined by 

Zr{G, q, v) = a-"(^)Z(G', q, v) = u^^^^Z{G, q, v) (2.25) 

which has the finite limit Zr — 1 as T — > 0. For a general strip graph {Gs)m of type Gg and 
length Lx = m, we can write 

ZriiG.U, q, a) = «^«^»)-) Yl ^o.,MG.,r = E cg.,,(Ag.,,,.)'" (2.26) 

with 

Ag.,,> = m^«^=)™)/"^Ag.j . (2.27) 

For example, for the strips of the square lattice with periodic longitudinal boundary 
conditions and free transverse boundary conditions, and of width Ly vertices, we have 



e{sq{Ly)m) = {2Ly — l)m, so the prefactor in ( p. 27] ) is 



u 



2L„-1 



III. ID CASE WITH FREE BOUNDARY CONDITIONS 



We first briefly discuss two cases that illustrate some important features in their simplest 
contexts. We begin with Potts/random cluster model on a line of n vertices, or, more 
generally, a tree graph, T„. One has the well-known result 

ZiT^,q,v) = qiq + vr-\ (3.1) 

This case illustrates the general feature that the antiferromagnetic random cluster model for 
real positive non-integral q fails to satisfy the usual statistical mechanical requirement that 
the partition function is positive, and hence does not, in general admit a Gibbs measure [p6| . 
Specifically, here we have 

Z(Tn,q,v)<0 for n even and q + v<0. (3.2) 

These negative values of Z{Tn,q,v) occur at physical finite temperature if < g < 1 and n 
is even, since the condition g + a — l<Ois equivalent to T < Tun, where 
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Tun = t-t4 7 = , < g < 1 , J < 0. (3.3) 

Although in this case one could restore positivity by requiring that n be odd, we shall show 
that there are further pathologies associated with this temperature. 

The continuous locus S = since the accumulation set of the zeros of Z in g is the 
discrete point q = —v. For q 0, the limits in the definitions (|2.11|) and ( p.lQ ) commute, 
and the free energy is (with v = a — 1) 

fgn = fnq = f = Hq + a-l) . (3.4) 

The physical thermodynamic behavior for this case will be compared below with that for 
the width Ly = 2 strips. For this purpose, we record the internal energy, 

Ja 

U = 3.5 

q+a-1 ^ ' 



and the specific heat. 



C^ '^^y-}^ . (3,6) 
(g + a - 1)^ 



Note that the specific heat ( p.6|) is positive if g > 1 but is negative and hence unphysical for 
all temperatures if g < 1, in both the ferromagnetic and antiferromagnetic cases. Thus the 
pathological nature of the range < g < 1 is manifested in the negative specific heat even 
for temperatures above the value Tun in eq. ( p.2| ) below which Z can be negative. Also, note 
that in the antiferromagnetic case there are unphysical divergences of U and C at T = T„„. 

For g = 0, the noncommutativity ( |2.8| ) occurs, and one has exp(/„q) = but exp(/g„) = 
(a — 1). This simple case demonstrates that the noncommutativity at a special point g^ can 
occur even when this point is not the the singular locus {B)qn or {B)nq- In passing we note 
that a study of the g-state Potts model on the Bethe lattice (tree graph in which the interior 



vertices all have the same coordination number) has been carried out in |57 . 



IV. ID CASE WITH PERIODIC BOUNDARY CONDITIONS 

A. General 

The Potts/random cluster model on the circuit graph C„, or equivalently, the ID line 
with periodic boundary conditions, is probably the simplest case with a nontrivial locus B. 
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The Tutte polynomial for this graph is well known |jT2|, and the corresponding Potts model 
partition function is 



Z(a,g,a) = (g + < + (g-l)^;". 



(4.1) 



As noted above, the Potts ferro magnet has a zero-temperature critical point, and this is 
also true of the antiferromagnet in the q = 2 case where these are equivalent. In the 
antiferromagnetic case, there is nonzero ground state entropy, S = kBlTi{q — 1) if g > 2. 

For comparison with the Ly = 2 results to be given below, we recall some of the thermo- 
dynamic properties of the ID Potts model. Here we take q > 2 where there is no pathological 
behavior (see below) and restrict to physical values of J and T; the resulting thermodynamic 
functions are then independent of whether one uses periodic or free boundary conditions and 
were given above for /, U, and C in eqs. (|3^ )- (|3.6| ). The internal energy and specific heat 
have the high-temperature expansions 



U = - 



J 



(4.2) 



(4.3) 



Recall that the T limit corresponds to K oo, i.e. u —>■ 0, for the ferromagnet (J > 0) 
and to K ^ — oo, i.e., a ^ 0, for the antiferromagnet (J < 0). The low-temperature 
expansions for these two cases are different: 



U 



-J 



1 - (g- l)e-^ +0(e 



~2K\ 



as 



K 



oo 



(4.4) 



U 



K r 



+ 0(e 



q-1 



2K\ 



as K 



-oo 



(4.5) 



C = kB{q - l)irV^ 1 - 2(g - l)e-^ + 0(e-'^) 



-K 



as K ^ oo 



(4.6) 



C 



q-1 



e-^ + 0{e 



2K\ 



as K ^ — oo . 



(4.7) 



Note that in the Ising case q = 2, these expansions satisfy the symmetry relations ( p.6| ) and 
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OTq ' ' i ' ' 2 ' ' 3 ' ' 4 ' ' 5 



K 



FIG. 1. Specific heat for the ID Potts ferromagnet as a function ot K = J/{kBT). Going from bottom 
to top in order of the heights of the maxima, the curves are for g = 2, 3, 4. 
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Olo~ ' ' i ' ' 2 ' ' 3 ' ' 4 ^ ^ 5 

-K 



FIG. 2. Specific heat for the ID Potts antiferromagnet as a function of —K = —J/(kBT). Going from 
top to bottom, the curves are for q = 2,3, 4. 

In Figs, m and |^ we plot C for the (ri — > oo limit of the) ferromagnetic (F) and anti- 
ferromagnetic (AF) cases (with fc^ = 1). In the antiferromagnetic case, C is a decreasing 
function of q for all < T < oo. In the ferromagnetic case, C increases (decreases) with q at 
low (high) temperatures and the curves for two different values q = qi and q = q2 are equal 
at the temperature T^r = J/ {ksKcr), where 

(4.8) 

Thus, Kcr — 0.35, 0.55, 0.90 for q = 2,3,4. The specific heat has a maximum at a temper- 
ature Kcmax given by the solution of the equation (g — 1){K + 2) + (2 — K)e^ = 0. Some 
illustrative values for Figs. |l| and |] are (i) Kcmax — ±2.4 (F,AF) for q = 2; (ii) Kcmax — 2.7 
(F), Kcmax ^ -2.2 (AF) for q = 3; (iii) Kcmax ^ 2.85 (F), Kcmax ^ -2.15 (AF) for q = 4. 
For q > 3, the value of C at this maximum increases (decreases) with q for the ferromagnet 
(antiferromagnet) . 



Kc 



A In 

2 



(gi-l)(g2-l) 
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B. l3g{{C}) for fixed a 



Returning to the study of Bq in the complex q plane as a function of a, we note that the 
solution of the degeneracy equation |g + a— 1| = |ci — 1| determines the locus Bq{{C}) to be 
the circle centered at g = 1 — a of radius |1 — a|: 

g = (l-a)(l + e^^) ,0<e<27r . (4.9) 

For a < 1 and a > 1, this circle has support in the right-hand and left-hand half planes 
Re{q) > and Re{q) < 0, respectively; for any a, it always passes through the origin. The 
locus Bg intersects the real q axis at g = and at g = qc{{C}), where 

q,{{C}) = 2(1 - a) . (4.10) 

In general, for finite n, the zeros of Z do not lie exactly on B. We have shown, however, that 
in the T = limit of the Potts antiferromagnet, i.e., for a = 0, the zeros of Z{Cn, q,a = 0) = 
P{Cn,q) do lie exactly on the locus B, which, for this case is the circle |g — 1| = 1 pTj^ . 
As T increases from to infinity for the Potts antiferromagnet, i.e. as a increases from to 
1, the radius and center of the circle both decrease from 1 to so that it contracts to the 
origin at a = 1. As a increases above 1 through real values, i.e. as T decreases from oo to 
for the ferromagnetic case, the circle is located in the Re{q) < half-plane, with its center 
moving leftward and its radius increasing as a function of a. In this ferromagnetic case, the 
crossing point given by qc{{C}) occurs on the negative real q axis; Bq does not cross the 
positive real q axis. 



C. BuiiC}) for fixed q 

We first consider values of g 7^ 0, 1, so that no noncommutativity occurs, and {Bu)ng = 
{Bu)qn = Bu- As discussed above, it is convenient to use the u plane since B^ is compact in 
this plane, except for the case g = 2, whereas Ba is noncompact because of the ferromagnetic 
zero-temperature critical point at u = 1/a = 0. For g 7^ 0, 1, 2, Bu is the circle 

Bu-. u = + e'^) , 0<cu<27r. (4.11) 

The exterior of this circle is the (complex-temperature extension of the) PM phase, and its 
interior is an O phase. 

For the ferromagnet, the fact that the singular locus Bu passes through the T = 
point M = for the ID Potts model with periodic boundary conditions, while for the same 
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model with free boundary conditions, Bu does not pass through m = 0, means that the 
use of periodic boundary conditions yields a singular locus that manifestly incorporates the 
zero-temperature critical point, while this is not manifest in Bu when calculated using free 
boundary conditions. As we shall show, this continues to be true concerning the longitudinal 
boundary conditions when one considers the Potts ferromagnet on the Ly = 2 strip graphs. 
This leads us to one of the important conclusions of this work, namely that although cal- 
culations of the free energy of the Potts model on infinite-length, finite-width strips with 
periodic boundary conditions in the longitudinal direction (the direction in which the strip 
length goes to infinity) are more difficult than with free longitudinal boundary conditions, 
the extra work is worth it since the resulting locus B incorporates this feature of passing 
through u = 0, corresponding to the zero-temperature critical point of the ferromagnet, if 
one uses periodic longitudinal boundary conditions. From our studies in the different, al- 
though related, context of chromatic polynomials [p^ , ^ , p] -|43[| , we reached the analogous 
conclusion that although the calculation of P{G, q) for strip graphs of a regular lattice is 
more complicated when one uses periodic longitudinal boundary conditions, the resultant 
singular locus Bq has the advantage of incorporating more of the features expected of the 
infinite- width limit, i.e. the full two-dimensional lattice. One such expectation, based on 
calculations of chromatic polynomials and the resultant W functions of eq. (|1 . 1 1|) for infinite- 
length, finite-width strips as the width increased, was that Bq passes through g = and that 
this nonanalytic locus separates the region including the interval < g < qc{{G}) on the 



real axis from the outlying region for sufficiently large |g| p2|j33|| ; this is also in agreement 
with the calculation in for the triangular lattice. However, for finite width strips, Bq 
consists of arcs ^2|] (and possible closed regions, as in Fig. 4 of |Q) which do not pass 
through g = and do not have this enclosure property. 
The circle ( [4.11|) crosses the real axis at m = and at 

•-«^» = MM) = -?^ 

(cf. eq. ( |4.1(]| )). The point Uc{{C}) occurs at complex temperature for q > 2 and physical 
temperature for q < 2. We shall comment further below on the case < g < 2. In the a 
plane, Ba is the vertical line 

Ba : Re{a) = Oc = 1 — ^, — oo < Im{a) < oo . (4-13) 

The phase with Re{a) > ac, to the right of this line, is the (complex-temperature extension 
of the) PM phase, while the phase to the left of the line is the O phase. As g — 2, the radius 
of the circle ( [4.11|) goes to infinity, and at g = 2, iS^ is identical to Ba by the symmetry 
relation ( p.l6| ), forming the full imaginary u axis. 
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We next consider the special values q = and 1 for which noncommutativity occurs. For 
g = 0, e^"'' = as in ( p.l3|) while in the PM phase with Re{a) > 1, e^^i" = a — 1 and in 



the O phase with Re{a) < 1, |e-^«"| = \a — 1|; the locus {Ba)nq = 0? while {Ba)qn is given by 
( [4.13| ) as the vertical line Re{a) = 1. For q = 1, the coefficient multiplying the second term 
in Z{Cn, q, v) vanishes, and Z{G, q = l,v) = a", a special case of ( p.2|) . Here e^"'' = a while 
in the PM phase defined by Re{a) > Oc = 1/2, we have e-^''" = a and in the O phase defined 
by Re{a) < 1/2, we have |e-^''"| = |a — 1|; {Ba)nq = since all of the zeros of Z occur at the 
discrete point a = 0, while {Ba)qn is given by eq. ( [4.13| ) as the vertical line Re{a) = 1/2. 

D. Phase Transition for Antiferromagnetic Case with < q < 2 

For the range < g < 2, and the antiferromagnetic case J < 0, the nonanalyticity in the 
free energy at a = = (2 — g)/2 in ([4.12|) occurs at the physical temperature 

Tp = 0<g<2. (4.14) 

Therefore, the generalization of the Potts antiferromagnet to real positive q defined by the 



random cluster representation (|1.8| ) has a finite-temperature phase transition in the n ^ oo 
limit of the circuit graph, i.e. in ID with periodic boundary conditions. (For the special 
value g = = 1, it is understood that one takes n oo ffist and then g — 1, i.e., 
one uses with the other order, g — > 1 and then — > oo, fnq is analytic and there is 
no phase transition.) The phase transition at T = Tp is not a counterexample to the usual 
theorem that a one-dimensional spin model with short-ranged interactions does not have any 
phase transition at finite temperature, because the existence of this transition is inextricably 
connected with the failure of positivity for Z and hence the absence of a Gibbs measure, 
which are implicit requirements for the applicability of the above-mentioned theorem. As g 
decreases from 2 to 0, the phase transition temperature Tp increases from to infinity. 

In the high-temperature paramagnetic phase T > Tp, the free energy, internal energy, and 
specific heat are given by the same expressions as for the n ^ oo limit of the Potts/random 
cluster model on the tree graph, eqs. (|3.4|), (P75|), and (|3.tiD, respectively. Hence, even in 
the high-temperature phase, one has an unphysical negative specific heat if g < 1. In the 
low-temperature O phase, strictly speaking, only \Z\ can be determined: | exp(/q„)| = |a — 1|, 
but with an appropriate choice of multiplicative phase, we can choose 

/ = ln(l - a) , T<Tp (4.15) 

and hence 
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U = , T<Tp (4.16) 

I — a 



C = -77^^ T<T,. (4.17) 



Thus, for all q in the range < q < 2 where there is a finite-temperature phase transition, 
the low-temperature phase has pathological property that the specific heat is negative. The 
phase transition itself is first-order, with latent heat 

lim U - lim U = . (4.18) 

T^T+ T^Tp q 

A basic pathology of the low-temperature phase of this antiferromagnet, i.e., the phase where 
|f I > |g -f- ^1, is that Z can be negative. For sufficiently large n, this occurs for 1 < g < 2 if 
n is odd and for < g < 1 if n is even. 

Thus, if one restricts to g > 1, this ID antiferromagnetic random cluster model satisfies, 
at least in the high-temperature phase, the requirement that the specific heat is positive and, 
for the interval 1 < g < 2, has a (first-order) finite-temperature phase transition; however, 
even if one restricts the approach to the n oo limit to even values of n, the low-temperature 
phase is unphysical because of the negative specific heat. One also observes that the results 
for the free energy and associated thermodynamic functions are the same for the n oo 
limit of the tree graph and the circuit graph, i.e. are independent of whether one uses free 
or periodic boundary conditions, if T > Tp, but differ for T < Tp, so that the existence 
of the low-temperature phase in the case of periodic boundary conditions also means that 
the n — > oo limit does not exist owing to different results obtained with different boundary 
conditions. The non-existence of a well-defined n oo limit for the random cluster model 



with non-integral g has been noted previously in ||56|. For positive integer g, the (zero-field) 
g-state Potts model is invariant under the operations of the permutation group Sg] however, 
this symmetry group is not defined for non-integral g. In any case, the usual Peierls argument 
shows that even if one could define some notion of a symmetry of Z for non-integral g, this 
symmetry could not be broken spontaneously in the phase transition for this ID system or, 
indeed, for the random cluster model on an infinite-length, finite-width strip, to be discussed 
below. 

A further generalization of this ID random cluster model is to keep g real but let it be 
negative; in this case the model with the ferromagnetic sign of the coupling, J > 0, formally 
has a nonanalyticity in the free energy at a positive finite value of the parameter T given by 
ksT = J/ln[(2 — g)/2]. However, since this model does not, in general, have a positive Z, 
one cannot really refer to this parameter as a physical temperature and we shall not discuss 
this case further. 
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E. Other slices of B{{C}) 



So far we have considered the "orthogonal" shoes of B obtained by holding either q or 
q constant. A different type of slice is obtained if one has q and a satisfy some functional 
relation. As an illustration of this, we consider perhaps the simplest case, namely the linear 
relation a + q = c, where c e C is a constant. If one treats the a and q variables as 
the "horizontal" and "vertical" axes (actually planes, in terms of real variables), then the 
condition a + = c is an affine translation of a diagonal slice of the complex locus B. The 
resultant Ba is the solution of the equation |c— l| = |a — 1|, which is a circle centered at 
a — 1 with radius |c — 1|. The corresponding Bq is a circle centered at q — c — 1 with radius 
|c — 1|. These circles in the a and q planes pass through a — and q — 0, respectively. 

V. SQUARE STRIP WITH FREE LONGITUDINAL BOUNDARY CONDITIONS 

In this section we present the Potts model partition function Z{Sm,q,v) for the Ly = 2 
strip of the square lattice Sm with arbitrary length = m + 1 (i.e., containing m + 1 
squares) and free transverse and longitudinal boundary conditions. One convenient way to 
express the results is in terms of a generating function: 



We have calculated this generating function using the deletion-contraction theorem for the 
corresponding Tutte polynomial T{Sm,x,y) and then expressing the result in terms of the 
variables q and v. We find 



oo 




(5.1) 



m=0 



r{S,q,v,z) 




(5.2) 



V{S,q,v,z) 



where 



X{S,q,v,z) 



As,o + As,iz 



(5.3) 



with 




(5.4) 




(5.5) 



and 



22 



V{S, q,v,z) = 1- {v-^ + Av^ + 3qv + q^)z +{v + l){v + qyv^z^ . (5.6) 

(The generating function for the Tutte polynomial T{Sm, x, y) is given in the appendix.) 
Writing 

2 

V{S,q,v,z) = \[{l-\s,,z) (5.7) 



we have 



where 



-^5,(1,2) — -j{Tsi2 ± \/ -R512 ) (5.^ 



^512 = v^ + Av^ + + g2 (5.9) 



and 



In we presented a formula to obtain the chromatic polynomial for a recursive family 
of graphs in the form ( |2.21| ) starting from the generating function. It will be useful to give 
here the generalization of this formula for the full Potts partition function. For a strip 
(recursive) graph with the labelling conventions used here, the generating function can be 
written as 

r(G, g, V, z) = — — (5.11) 



with 



V{G,q,v,z) 



U{G,q,v,z) = Y,^G,jZ^ (5.12) 
i=o 



and 



V{G,q,v,z) = l + J2bGj 



■z^ 



d-D 



11(1 -Xa,z) (5.13) 



where 
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d^f{G) = deg,{M{G)) 



(5.14) 



dv{G) = deg-XV{G)) 



(5.15) 



Then the formula is 



Z{Gra,q,v) 



s=0 



n 



1 



\ m 



(5.16) 



For our present open strip Sm, we have 

(^5,0-^5,1 + As,l] 



Z{Sm, q,v) 



(-^5,1 — -^5,2) 



A 



5,1 



_l_ (^5,oAs,2 + ^5,1) ^ 
(A5,2 — A^^i) 



S,2 



(5.17) 



(which is symmetric under Xs^i ^ ^8,2)- This shows that the cgj can depend on both q and 
V for open strip graphs. Although both the A5j's and the coefficient functions involve the 
square root y/Rsi2 and are not polynomials in q and v, the theorem on symmetric functions 
of the roots of an algebraic equation [^] guarantees that Z{Sm, q,v) is a polynomial in q 
and V (as it must be by ( |1.8| ) since the coefficients of the powers of z in the equation ( p. 131 ) 
defining these A^^'s are polynomials in these variables q and v. This is a generalization of 
our discussion in |^1[] from the special case of chromatic polynomials to the general case of 
the Potts/random cluster partition function. 

As will be shown below, the singular locus Bu consists of arcs that do not separate the 
u plane into different regions, so that the PM phase and its complex-temperature extension 
occupy all of this plane, except for these arcs. For physical temperature and positive integer 
q, the (reduced) free energy of the Potts model in the limit n — ^ 00 is given by 



/ 



InA 



{5.U 



This is analytic for all finite temperature, for both the ferromagnetic and antiferromagnetic 
sign of the spin-spin coupling J. The internal energy and specific heat can be calculated 
in a straightforward manner from the free energy ( p. 18] ); since the resultant expressions are 
somewhat cumbersome, we do not list them here. We find that for g < 2, in both the 
ferromagnetic and antiferromagnetic case, for sufficiently low temperature, the specific heat 
is negative, and hence the random cluster model is unphysical for g < 2 on this family of 
graphs. 

Let us define 



Mq) 



P{Ck,q) 

q{q - 1) 



k-2 



-11 



s 



k-2-s 



(5.19) 
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and P{Ck, q) is the chromatic polynomial for the circuit (cyclic) graph Ck with k vertices, 

P{Ck^q) = {q-lf + {q-l){-lf (5.20) 
so that D2 = I, = q — 2, 

Di = q'^-3q + 3 (5.21) 

and so forth for other D^s. In the T = Potts antiferromagnet limit v = —1, A^-^i = D4 and 
^s,2 = 0, so that eq. ( |5.2|) reduces to the generating function for the chromatic polynomial 
for this open square strip (cf. eq. (2.16) in p^ ) 



T{S,q,v = -l;z) = — 5.22 

1 — D4Z 

where Equivalently, the chromatic polynomial is 

P{Sm,q)=qiq-l){D,r^' . (5.23) 

For the ferromagnetic case with general q, in the low-temperature limit v 00, 

= t;3 + 3^2^ (^^^2)^7 + 0(1) , \s,2 = v^ + 2{q-l)v + 0{l) as v ^ 00 (5.24) 

so that \Xs^i\ is never equal to |A5^2| in this limit, and hence Bu does not pass through the 
origin of the u plane for the n ^ 00 limit of the open square strip: 

u = O^Bu{{S}). (5.25) 

In contrast, as will be shown below, Bu does pass through -u = for this strip with cyclic 
or Mobius boundary conditions. For our later discussion, we record here the expressions for 
the Xs/s for the Ising case, q = 2: 

A5,(i,2) = liv + 2)\v^ + 2v + 2± {v^ + Av^ + 8w + 4)^/^1 _ (5_26) 



A. Bq{{S]) for fixed a 

We discuss here the continuous locus Bq{{S}) in the q plane for various values of a. For 
the chromatic polynomial case a = (w = —1), S = 0, since W{{S},q) = (1^4)^/^ has only 
the discrete branch point singularities (zeros) at 

q^p, g^p = l + e±*-/3_ (5.27) 
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However, for a 7^ 0, the situation is qualitatively different; Bq is nontrivial. As a increases 
above 0, the locus Bq forms two complex-conjugate (c.c.) circs, cLS shown, for a = 0.1, in Fig. 
^. For small a, these arcs lie near to the circle ji? — 1| = 1; as a decreases, they shorten and as 
a — > 0, they degenerate to the c.c. points qbp, ql^ in eq. (|5.27| ). The endpoints of the arcs are 
the (finite) branch point singularities of \s,ji j = 1;2, arising from the zeros of the square 
root in (|5.8| ); for example, for a = 0.1, these endpoints occur at g ~ 1.0654 + 0.9293i and 
q ~ 1.6346 + 0.59275i, together with their complex conjugates. From Fig. ^ one can see that 
the density of chromatic zeros is greatest at the endpoints and minimal at the centers of the 
arcs. As a increases, these arcs extend downward toward the positive real q axis. As a reaches 
the value a = 9/16, the arcs touch the real axis at g = 63/64 = 0.984375, thereby joining to 
form a single self-conjugate arc with endpoints at q,q* = (21± 14v/6i)/64 ~ 0.3281±0.5358i 
As a increases above the value 9/16, Bq consists of the self-conjugate arc and a line segment 
on the real axis, which spreads out from the point q = 63/64. This corresponds to the fact 
that for a > 9/16, the expression in the square root in eq. (|5.^ ) has real as well as complex 
zeros. An illustration is given for a = 0.9 in Fig. ^. As a ^ 1, the locus Bq shrinks in toward 
the origin. For the ferromagnetic range a > 1, i3g is located in the Re{q) < half-plane and 
forms c.c. arcs together with a line segment on the negative real illustrated for a = 2 

in Fig. |[ 
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FIG. 3. Zeros of Z{Srmq,n) in the q plane for a = 0.1. For this and the other figures on zeros of 
Z{Sm, q, a), we use m = 19, i.e., n = 42. The axis labels are qr = Re{q) and qi = Im{q) here and in other 
g-plane plots. 
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FIG. 4. Same as Fig. | for a = 0.9. 
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FIG. 5. 


Same as Fig. S for a — 2. 





One can also consider negative real values of a, which correspond to complex temperature. 
As a decreases from through real values, Bq forms arcs, as was the case when a increased 
from 0; these arcs have endpoints at the branch point singularities of Xsj and elongate as a 
moves downward in the range — l<a<0. As a reaches —1, these arcs touch the positive 
real q axis at the point q = 2 and join to form a single self-conjugate arc (with endpoints 
at 4 ± 4z), but as a decreases below —1, the arcs retract from the real axis to form two c.c. 
parts again. It is straightforward to consider complex values of a also, but we shall restrict 
ourselves to real a here. 

B. Bu{{S}) for fixed q 

For our analysis of Bu{{S}) we start with large q. Here Bu consists of a self-conjugate arc 
that crosses the real u axis, together with a complex-conjugate pair of arcs that are concave 
toward the real u axis. As g ^ oo, these arcs all shrink and move in toward the origin of 
the u plane. This limit thus commutes with the result of taking g ^ oo first before taking 
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n — >■ oo; in this case, 

As,i = + 3wg + 0(1) , Xs,2 = v\l + v) + O (^^^ as q oo (5.28) 

so that the degeneracy equation jA^^il = |A5'^2| has no solution for g — oo. In Fig. || 
we show the complex-temperature zeros of Z for a typical value, q = 10, calculated for 
Lx = m + 1 = 20, i.e., n = 42. With this large a value of n, these zeros occur close to the 
asymptotic locus Bu and give an adequate indication of its location. The self-conjugate arc 
crosses the real u axis at m ~ —0.3954 where the quantity T512 in eq. (|5.9|) vanishes, so that 
|As,i| = 1^5,21- The endpoints of this arc occur at two of the zeros of the square root in eq. 
i f).8\ ), viz., u ~ — 0.2937 ± 0.3870i. The two c.c. arcs have their endpoints at the four other 
zeros of this square root, at m ^ —0.1361 ± 0.14245i and u ~ 0.1178 ± 0.8130i. 

As q decreases, the endpoints of the self-conjugate arc retract toward the real axis, it curls 
over to be more concave to the right, and the point at which it crosses the real axis moves 
to the left. For example, for g = 4, the self-conjugate arc crosses the real axis ai u = —1 and 
has its endpoints at g ~ —0.4341 ± 1.3178i, and the c.c arcs extend between endpoints at 
q ~ -0.3697 ± 0.2394Z and q ^ 0.1711 ± 0.1593i. For g = 3 (see Fig. 0), the self-conjugate 
arc crosses the real axis at u ~ —1.2767 and has endpoints at g ^ —0.5498 ± 0.2489i, and 
the c.c. arcs extend between g ~ —0.0839 ± 2.0177^ and 0.1892 ± 0.1974i, passing through 
the points u,u* = e^^'^*/^. These results for g = 3 and 4 have interesting implications that 
we shall discuss further below. The changes in Bu as g decreases further toward g = 2 are 
illustrated in Fig. |] where we show the g = 2.5 case. The self-conjugate arc crosses the real 
axis at M ~ —1.323 and has endpoints at m ~ —0.7249 ± 0.2083^ while the c.c pairs of arcs 
have endpoints at g ^ 0.2006 ± 0.2272z and g ~ 0.56505 ± 2.4351i. 

For g = 2 (Fig. |^), the quartic polynomial in the square root of eq. ( p.8| ) factorizes 
into a quadratic polynomial times {v + 2)^, yielding the result (|5.26|) . Correspondingly, 
the self-conjugate arc disappears, and the locus B consists of two complex-conjugate arcs 
located in the half-plane Re{q) > 0, and touching the imaginary axis at m = ±i. This locus is 
invariant under the inversion map u ^ 1/u. The upper arc extends from the left endpoint at 
Mel ^ 0.2138 + 0.27202 through u = i to a right endpoint at Ue2 = 1/m*i ^ 1.78615 - 2.2720z, 
and so forth for the c.c. arc. As we have discussed before in the context of Bq (see 
also |]59|j60[| ), these endpoints are the zeros of the square root in (|5.26| ) where there are 
finite branch point singularities in A^^i. There is also a discrete zero of Z at the point 
u = —1 with multiplicity scaling proportional to the lattice size. As we proved in a previous 
theorem (Theorem 6 of Ref. |^), this zero arises for the Ising model on a lattice with odd 
coordination number; in the present case, all of the vertices of the strip S except those on 
the four end-corners have A = 3. 
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Because of the unphysical nature of the Potts/random cluster model for g < 2, we shall 
not discuss this range except to mention another example of the noncommutativity (|2.8| ) at 
q = 0. If one first sets q = and calculates Z, then, since Z = identically, the set of zeros 
of Z is vacuous. However, if one takes the limit n — oo, calculates the accumulation set Bu, 
and then takes the limit g ^ 0, one finds that limg^i Bu is not the empty set. This is clear 
from the fact that in the limit g — 



A 



S,{l,2);q=0 



(5.29) 



so that the degeneracy equation |A5,i;g=o| = |A5',2;g=o| has a nontrivial solution, namely the 
section of the the circular arc in the Re{u) < half-plane 

Bu-. u = ie*^''^^) , 0<9< arctan{2V2) (5.30) 
o 

which crosses the real axis aX u = —1/3 (i.e., v = —4) and has endpoints at u = (—1 ± 
2\/2i)/9. 



31 



. 6-- 




-0.6- 

-1 ' -o'.'s' -o'.'e' '-'o'.V '-'o'.2' ' 'i' ' ' bli' 'oik' 'ole' '6.8 

ur 



FIG. 6. Zeros of Z{Sm,Q,0') in the u plane for g = 10 and m = 19 (n = 42). The axis labels are 
ur = Re{u) and ui = Im{u) here and in other M-plane plots. 
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FIG. 7. Same as Fig. || for g = 3. 
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FIG. 8. Same as Fig. |6|for q = 2.5. 
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FIG. 9. Same as Fig. || for g = 2. 

VI. CYCLIC AND MOBIUS LADDER GRAPHS 
A. Results for Z 

By either using an iterative application of the deletion-contraction theorem for Tutte 
polynomials and converting the result to Z, or by using a transfer matrix method (in which 
one starts with a. x transfer matrix and generalizes to arbitrary g), one can calculate 
the partition function for the cyclic and Mobius ladder graphs of arbitrary length, Z{G, q, v), 
G = Lm, MLm- We have used both methods as checks on the calculation. Our results have 
the general form ( [^.18D with Nx = 6 and are 



Z(L^,g,t;) = 5^c^,,(Az.,,r (6-1) 



and 
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where 



Z{MLm,q,v) = ^CML,j{>^ML,j 



^ML,j — ^L,j , j — 1, •••5 6 



A. 1 = 



Al,2 = v{v + q) 



A 



i,(3,4) 



V 

2 L 



q + v{v + 4:)± (f ^ + 4f ^ + 12t;^ - 2gt;^ + Aqv + g' 



2N1/2 



and 



A 



L,5 



As,! ; A 



L,6 — Ag 2 



(6.2) 

(6.3) 
(6.4) 
(6.5) 
(6.6) 

(6.7) 



where A^j, j = 1,2 for the open ladder were given above in eq. (p.8|). We note that 
Al,3Al,4 = {l + v){q + v)v^ and Xl,5^l,6 = (1 + f )(q' + f )^f ^. Chromatic and Tutte polynomials 
for recursive families of graphs obey certain recursion relations p2| , |3^ . In terms of the 
equivalent Tutte polynomial, given in the appendix, the results ( |6.1| ) and ( |6.2| ) agree with a 
recursion relation given in Ref. |^2|] (see also p3|). 

The coefficient functions for the cyclic and Mobius ladders are 



cl,i = - 3g + 1 



Cl,2 — Cl,3 — Cl,4 — CmL,3 — CmL,A — 9 " 1 



Cl,5 — CLfi — CmL,5 — CmL,6 — 1 



cmla 



Cml,2 = 1 - q . 



(6.8) 
(6.9) 
(6.10) 
(6.11) 
(6.12) 



Because of the equalities cg,3 = cg,4 and c^.s = cg,6 for G = L and for G = ML, we can 
again apply the theorem on symmetric polynomial functions of roots of algebraic equations 
|6l| to confirm that, although the Ag-j's for nonpolynomial algebraic functions of q and v for 



j = 3,4,5,6, Z{Gm,q,v) is a polynomial function of these variables q and v, as it must be 
by O- 
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B. Special values and expansions of A's 



We discuss some special cases. First, for the zero-temperature Potts antiferromagnet, i.e. 
the case a = {v = —1), the partition functions Z{Lm,q,v) and Z{MLm,q,v) reduce, in 



accordance with the general result (|1.1CI|) , to the respective chromatic polynomials P{Lm, q) 
and P{MLm,q) calculated in ||22|. In this special case, we have Al^i = 1, Xl,2 = ^ — q, and 
(for an appropriate choice of sign of terms of the form \/{q — 3)^ and ■>/ (-04)^ ) 3 = 3 — g, 
Al,4 = 0, Xl,5 = -D4 = — 3g + 3, and = 0. For the infinite-temperature value a = 1, 



we have Xlj = for j = 1,2,3,4,6, while Xl,5 = q"^, so that Z{G,q,a = 1) 
G = Lm, MLm, in accord with the general result ( p.3|) . 
At g = 0, besides the g-independent A^^i, we find 



q 



2m 



A 



L,2 



g" for 



(6.13) 



A 



L,3 



A 



L,5 



(a -11 



a + 3 + (a^ + 2a + 9) 



1/2 



(6.14) 



-(a - I] 
2^ 



a + 3 - (a^ + 2a + 9) 



1/2 



Al,4 — Xlfi 

Since A^^s and A^^s are leading and are degenerate at this point, it follows that 

g = is on Bq{{L}) V a . 



(6.15) 



(6.16) 



At g = 1, CL,j = for j = 2,3,4 so that the corresponding Xlj, j = 2,3,4, do not 



contribute to Z. Further, cl,i 



-CL,j, j = 5,6 and A^^i 



A, 



so 



that the contributions of these terms cancel in Z, leaving only the contribution of A^ 5: 



Z{Lm,q = I, a) = (Al,i 



,3m 



in agreement with the general formula 



In order to study the zero-temperature critical point in the ferromagnetic case and also the 
properties of the complex-temperature phase diagram, we calculate the Acj.m's corresponding 
to the Agj's, using eq. ( ^^.271 ). This gives Xl,i,u = m(1 - m)^, Xl,2,u = u{l - ■u)[l + (g - 1)m)], 
and so forth for the others. In the vicinity of the point u = one has 



Xliu = u — 2u + u 



(6.17) 



AL,2,n = M + (g - 2)u^ + (1 - q)u^ 
and the Taylor series expansions 



(6.18) 
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Al,3,« = 1 - + 2(g - 2y + 0{u^) (6.19) 
Al,4,« = u+{q- 4)u^ + (7 - 3g)M3 + O(m^) (6.20) 
AL,5,n = 1 + (g - l)w' [l + 4m + 0(^2)] (6.21) 

Al,6,« = m + 2(g - 2)^2 + (g2 _ 7g + 7)^3 ^ ^^^4^, _ 22) 

Hence, at m = 0, Xl,3,u and A^^s^^ are dominant and |Al, 3,^(1 = |Al,5,u|, so that the point u = 
is on Bu for any q 7^ 0, 1, where the noncommutativity ( |2.8| ) occurs. For g > 0, Al^5_„ is 
dominant on the real u axis in the vicinity of m = and hence in the PM and O phases that 
can be reached by analytic continuation therefrom, while the term XL,3,r is dominant on the 
imaginary u axis in the neighborhood of the origin, and hence in the O phases that can be 
reached by analytic continuation from this neighborhood. 

To determine the angles at which the branches of Bu cross each other at m = 0, we write u 
in polar coordinates as m = re*^, expand the degeneracy equation |A2,^3^„| = |Al^5_„|, for small 
r, and obtain gr^ cos(2^) = 0, which implies that (for g 7^ 0, 1) in the limit as r = \u\ 0, 

e=^ll±}^, J = 0,1,2,3 (6.23) 

or equivalently, 6 = ±7r/4 and 6 = ±37r/4. Hence there are four branches of Bu intersecting 
at M = and these branches cross at right angles. The point m = is thus a multiple point 
on the algebraic curve Bu, in the technical terminology of algebraic geometry (i.e., a point 
where several branches of an algebraic curve cross ||62|| ). 

In order to investigate how these crossings depend on Ly, we have calculated Z for the 
cyclic strip graph of the square lattice with the next larger width, Ly = 3. Since the T = 
critical point for the Potts ferromagnet is present for each q 7^ 0, 1, it suffices to do this 
calculation for the simple q = 2 Ising case (bearing in mind the noncommutativity that 
applies at special values as discussed above). We find that there are two Aj's that are 
dominant near u = 0, and the small-u expansion of the degeneracy equation yields the 
condition cos(3^) = 0, so that there are six curves on Bu crossing m = 0, at the angles 
6 = {2j + l)vr/6, j = 0,1,. ..,5. This leads to the generalization that for the cyclic strip 
graph of the square lattice with width Ly, there are 2Ly curves on Bu that cross each other 
at M = 0, at the angles 9 = (2j + l)n/{2Ly), j = 0, 1, ...,2Ly — 1. This inference implies, 
in turn, that in the limit Ly —>■ 00, an infinite number of curves on Bu intersect at u = 0, 
and the complex-temperature (Fisher) zeros become dense in the neighborhood of this point. 



38 



Since the origin of this phenomenon is not dependent in detail on the lattice type, one would 
also infer that it occurs for infinite-length width Ly cyclic strips of other lattices. 

For g = 0, 1 the ( pl2D occurs, with {Bu)nq = by eqs. (|23|) and (|2J[5| ), but {8u)qn ^ 0- 



While Bu is compact for g 7^ 2, it is noncompact for q = 2, where the symmetry (|2.16|) holds. 
Our exact calculations yield the following general result 

B{{L}) = B{{ML}) . (6.24) 

This is in accord with the conclusion that the singular locus is the same for an infinite- 
length finite-width strip graph for given transverse boundary conditions, independent of 
the longitudinal boundary condition. This generalizes our previous finding that Bq was 
independent of the longitudinal boundary conditions for the case a = p7| , p| , ^ . In the 
present case, the result ( |6.24| ) follows immediately because Z{Lm,q,v) and Z{MLm,q,v) 
involve the same A^'s. We note that this is a sufficient, but not necessary condition for 
the loci to be the same for a given family of graphs when one changes the longitudinal 
boundary conditions; it may be recalled that for the T = Potts antiferromagnet on the 
width Ly = 3 strip of the square |^ with periodic transverse boundary conditions, when 
one changed from periodic to twisted periodic longitudinal boundary conditions, i.e. toroidal 
to Klein bottle topology, three of the A^a = 8 terms were absent. However, since none of 
these was a dominant term anywhere, the locus Bq was the same for either toroidal or 
Klein bottle boundary conditions. From our calculation of the chromatic polynomial for the 
width Ly = 3 strip of the triangular lattice with both free and periodic transverse boundary 
conditions and periodic and twisted periodic longitudinal boundary conditions |^ we found 



that a similar situation occured for the toroidal versus Klein bottle boundary conditions: six 
of the A^A = 11 terms in the toroidal case were absent in the Klein bottle case, but again 
none of these was dominant anywhere. Owing to the equality (|6.24]) , we shall henceforth, for 
brevity of notation, refer to both B{{L}) and B{{ML}) as B{{L}) and similarly for specific 
points on B, such as qc{{L}) = qc{{ML}), etc. 

C. Bq{{L}) for fixed a 

We find that Bq{{L}) crosses the real q axis at 

= (l-a)(a + 2) . (6.25) 

This is the solution to the degeneracy equation of leading terms lA^^sl = \\l,3\ = |'^l,2|- As a 
increases from to 1, gc({-^}) decreases monotonically from 2 to 0. From eq. (|6.25|) it follows 
that there are, in general, two values of a that correspond to this value of q on B{{L}), viz.. 
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«c,±(W) = ^[-l±v/9^], ^.e, u^A{L})= ^^(g^y^ • (6-26) 



D. Antiferromagnetic Case, T = 

We start with the T = antiferro magnet, i.e. the case a = 0. After initial studies in 
Refs. P^J2B|J25[ ], the locus B was determined in Ref. [^]. As is shown in Fig. 3 of Ref. 



Eq separates the q plane into four regions, and qc{{L}) = 2. The outermost region is -Ri, and 
includes the segments 2 < q and g < on the real q axis; in this region A^^s is dominant. 
The innermost region, denoted R^, includes the segment < g < 2 on the real axis; in 
this region, the term Al^3 is dominant. In addition, there are two other complex-conjugate 
regions, R2 and R2, which touch the real axis at g = qc{{L}) = 2 and stretch outward to 
triple points at 

QLMP., Ct-P. = 2±v^^- (6.27) 



The part of B separating region R^ from regions R2, R2 is the line segment Re{q) = 2, 
-V2 < Im{q) < V2. In regions R2, R2, ^2 is dominant. At g = 2 all four terms are 
degenerate (recall that for a = 0, \l,4 = Xl,6 = 0). At the triple points qL,trip., there are 
three degenerate leading terms, with lA^^sl = |Al,3| = |Al,2|- All four regions are contiguous 
at qc{{L}). 

E. Antiferromagnet Case for T > 

We proceed to consider the regions in the g plane for the Potts antiferromagnet at ar- 
bitrary nonzero temperature, i.e. the range < a < 1. The zeros of Z in the g plane are 
shown for several values of a in the figures. In this range < a < 1 we find a number of 
general features. As was true at T = 0, Bq continues to separate the g plane into different 
regions and, as indicated in eq. ( |6.16| ) and ( p.25| ), this locus crosses the real axis at g = 



and qc{{L}). Bq consists of a single connected component made up of several curves. Com- 
menting on the regions in the g plane, starting for a near 0, we note that again the region 
Ri is the outermost, and includes the semi-infinite line segment on the real axis g > qc{{L}) 
and g < 0; region R3 is the innermost region, and includes the line segment < g < qc{{L}). 
The complex-conjugate regions R2 and R2 extend upward and downward from gc({L}) to 
triple points. As a increases, the complex-conjugate regions R2 and R2 are reduced in size. 
As is evident in the figures, as a increases from to 1, the locus Bq contracts toward the 
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origin, q = and in the limit as a — 1, it degenerates to a point at g = 0. This also describes 
the general behavior of the partition function zeros themselves. That is, for finite graphs, 
there are no isolated partition function zeros whose moduli remains large as a — 1. This is 
clear from continuity arguments in this limit, given eq. ( p.3|) . 




FIG. 10. Zeros of Z{Lm, q, a) in the q plane for a ~ 0.25 and to = 18 (n = 36). 
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FIG. 11. Same as Fig. ITolfor a = 0.5 
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FIG. 12. Same as Fig. IT^ for a = 0.75. 



F. Ferromagnetic Case 



In Fig. |1^ we show the zeros of Z for a typical ferromagnetic values, a = 2. We find 
the following general features of Eg for the ferromagnetic Potts model for the full range of 
temperature. The locus Bg contains a heart-shaped figure and a finite line segment on the 
negative real q axis. The line segment occurs because the expression in the square roots 
in Al,5 and Xifi, given as Rsi2 in eq. (|5.10| ), is negative in an interval of the negative real 
axis, yielding a pure imaginary square root so that, given that Tsu is real, |Al,5| = |Al,6|- 
For example, for the case shown in Fig. Rsi2 < for —3.73 < q < —2.10; within this 
interval, lA^^sl = |Al,6| are leading for —3.73 < q < 3.35, thereby producing the line segment. 
(For the remaining part of the interval, —3.35 < q < —2.10, these eigenvalues have smaller 
magnitudes than lA^al and hence do not determine the locus Bg.) The size of the heart- 



shaped boundary increases as a increases. Since qc{{L}), given in eq. (|6.25| ), is negative, 
Bg does not intersect the positive real q axis. As was true for the Potts AF, in the region 
exterior to Bg in the q plane, the dominant \j is A^^s, so that the (reduced) free energy is 
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(6.28) 



where Xs^i was given in eq. (^). For the range q > 2 where this system has acceptable 
physical behavior, the above expression for the free energy holds for all (physical) tempera- 
tures. We shall discuss the thermodynamics further below. In the region interior to Bq, Xl ^ 
is dominant, so |e-^| = lA^^sl"^/^. 




FIG. 13. Zeros of Z{Lm, q, a) in q for a = 2 and m = 18 (n = 36). 



G. Bg for a < 

We briefly comment on Bq for negative real values of a, which correspond to complex 
temperature (as well as complex values of a not considered here). For the interval 

1 

- - < a < (6.29) 

the regions R2 and R2, which were separate, although contiguous at q = qc{{L}), for a in 
the interval < a < 1, now merge to form one region, which we shall call R22* to indicate 
this merger. In this region, 2 is dominant. The R22* — Ri boundary is determined by 
the degeneracy equation |Al,2| = |Al,5| and crosses the real axis at qc{{L}). The point at 
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which the — R22* boundary crosses the real axis is determined by the relevant root of the 
threefold degeneracy equation \\l,3\ = \^l,2\ = I-^l.sI; and is 



Qbi = 2(1 — a) for — - < a < . 



(6.30) 



The width of the merged R22* region on the real axis is thus (— a)(l — a) for this range of 
a. As a decreases in the interval ( |6.29|) , qhi and qc both increase above 2. When a decreases 
through the value a = —1/2, at which point qbe = 9/4, the square root in X^ ^ becomes 
complex. Viewed the other way, solving eq. ( |b.3(JD for a gives 



_ 1 



1 ± V9 - 4g 



(6.31) 



We are interested in the larger solution, ab+. When q increases through 9/4, corresponding 
to ab+ = —1/2, the square root becomes complex, and there is no longer a real solution for 
ab±. The region diagram changes qualitatively for a < —1/2. The illustrative case a = —1 
is shown in Fig. |14|. Note that Z{Lm,q,CL = —1) has an overall factor q{q — 2). Here Bg 
crosses the real q axis at g = 2 and g = gc = 4 as well as at g = 0. The crossing at g = 4 is 
a multiple point on the algebraic curve. Other interesting changes occur for larger negative 
values of a, but we shall forgo discussing them to proceed to the physical range of real a > 1. 




FIG. 14. Zeros of Z{Lm, q, a) in g for a = —1 and m — 18. 
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H. Thermodynamics of the Potts Model on the Ly = 2 Strip 

1. Q>2 



In this section we first restrict to the range q >2 where the Potts/random cluster model 
has physical behavior for both the ferromagnetic and antiferromagnetic cases, and then 
consider the behavior for < g < 2. For q > 2, the free energy is given for all temperatures 
by ( |6.28| ). It is straightforward to obtain the internal energy U and specific heat C from this 
free energy; since the expressions are somewhat complicated we do not list them here but 
instead concentrate on their high- and low-temperature expansions and general features, as 
compared with those for the Ly = 1 case. The high-temperature expansion of U is 



U 



3J 
2^ 



q 



(6.32) 



The expression in brackets is the same as that for the Ly = 1 strip up to and including the 
term. For the specific heat we have 



C 



3k 



2q- 



(6.33) 



(here the order term differs from that for Ly = 1.) The low-temperature expansions for 
the ferromagnet {K oo) and antiferromagnet {K — oo) are 



U = J 



-- + (g - l)e-^^ 1 + 6e-^ + 7{q - l)e-'^ + 0{e 



and 



U 



2{D,y 



h + t2e^ + 0(e 



-,2K\ 



as K 



-oo 



as K ^ oo (6.34) 



(6.35) 



where 1^4 = — 3g + 3 was given in ( ^.21| ) and 

ti = (g-2)(3g2-9g + 8) 



(6.36) 



(3g6 - 42g5 + 211g4 - 532g3 + 734^2 _ 534^ + 162) 

GDI? 



(6.37) 



C = 2kBK\q - l)e-'^ 1 + 96"^ + 14(g - 1)6"'^ + 0(e 



as K 00 



(6.38) 



and 



46 



c 



ksK'e 



as 



K 



-oo 



(6.39) 



Again, we observe that for the Ising case q = 2, these expansions satisfy the symmetry 
relations ( |2.6| ) and ( |2.7| ). (In passing, we mention the generahzation of the first term in eq. 
( |S.34| ) to arbitrary Ly-. in the T = hmit of the Potts ferromagnet. 



U 



2Ly\ 



J . 



(6.40) 



We show plots of C (with fc^ = 1) for the ferromagnetic and antiferromagnetic Potts model 



on the Ly = 2 strip (in the oo limit) in Figs. O and ITBl As was true for Ly = 1, in 



the antiferromagnetic case, C is a decreasing function of q for all finite temperature, while 
in the ferromagnetic case, C increases (decreases) with q at low (high) temperatures and 
has a maximum that increases with q. For a fixed q, by comparing the previous plots of the 
specific heat on the line {Ly = 1 case) with the corresponding plots for the Ly = 2 strip, for 
the ferromagnet, and for the antiferromagnet, one can see quantitatively how the behavior 
of this function changes as Ly increases. 




FIG. 15. Specific heat for the Potts ferromagnet on the infinite-length, width Ly ^ 2 strip (ladder) as a 
function of K — J/{kBT). Going from bottom to top in order of the heights of the maxima, the curves are 
for q = 2,3,4. 
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-K 



FIG. 16. Specific heat for the Potts antiferromagnet on the infinite-length, width Ly = 2 strip (ladder) 
as a function of —K = —J/{kBT). Going downward in order of the heights of the maxima, the curves are 
for g = 2, 3, 4. 

For both the Ly = 1 and Ly = 2 strips, we observe that the exponential zero in the 

specific heat as T — > for both the ferromagnetic and antiferromagnetic cases is C ~ 

{q — l)i<'^e~^^l^l. For comparative pm poses we have also calculated the partition function, 

free energy, and these thermodynamic functions for the Ising model on the strip with the 

next larger width, Ly — 3. We find that the above dependence on Ly is again exhibited, 
namely C ~ X2g-3|Jf I 

In view of the fact that the Potts ferromagnet has a zero-temperature critical point 
for the infinite-length, finite-width strip graphs of the square lattice (as does the Potts 
antiferromagnet in the q — 2 case where it is equivalent to the ferromagnet on these graphs) , 
it is of interest to investigate the dependence of the singularities in thermodynamic functions 
on the strip width Ly. As is typical for systems at their lower critical dimensionality, these 
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are essential singularities. We have done this comparative study above for the specific heat. 
We next consider the (divergent) exponential singularities in the correlation length and the 
(uniform, zero-field) susceptibility (again, in the q = 2 case, under the replacement J ^ — J 
and uniform staggered, this subsumes the antiferromagnetic case). Define the ratio of the 
subleading eigenvalue divided by the leading eigenvalue of the transfer matrix for the strip 
graphs with periodic longitudinal boundary conditions considered here as piy- Thus 

Pi = - — = —— 6.41 

Xc,i q + v 

and 

P2 = ^ . (6.42) 

We have also calculated ps for the Ising case, but since it is a rather messy expression 
involving cube roots, we do not display it here. These ratios control the asymptotic 
decay of the spin-spin correlation function. For example, in the n —>■ oo limit, the spin-spin 
correlation function in the ID case is given by {arCr') oc (pi)!*""^ '. The correlation length 
can be written as 

e = . (6.43) 

In PL, 

For the ID case, one knows that the correlation length has an exponential divergence as 
T ^ 0: ^ ~ q-^e^ + 0(1). For the Ly = 2 strips we find 

^ ~ g^^e^^' + O(e^) , as T ^ for {G} = {L} . (6.44) 

and for the Ising model on the width Ly = 3 strip, we obtain ^ ~ (l/2)e^^ + 0(e^'^). These 
results show that the exponential divergence in the correlation length is more rapid for larger 
width Ly and are consistent with an inference that ^ ~ q~^e^y^ + 0{e^^y~^^^) as T — 0. 
The fact that the correlation length diverges more rapidly as Ly increases is easily explained 
since this is due to the spin-spin interactions and the average effect of these interactions, as 
determined by the average coordination number, Aave in eq. ( p.5|) , increases as Ly increases. 

The zero-field susceptibility (per site) is well known for the ID case: x = Pi^+Pi) I (l~Pi)) 
which diverges as a function of K like x ~ Ke^ as K ^ oo. Our results for Lj^ = 2, 3 support 
the inference that x ~ Ke^y^ as K ^ oo. The more rapid divergence in x as Ly increases 
can be explained in the same way as was done for the correlation length. 

The inferred Ly dependence of the divergences in the correlation length ^ and susceptibil- 
ity X a-t the zero-temperature critical point of the Potts ferromagnet dramatically illustrate 
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the fact that the thermodynamic behavior of the model on this sequence of infinite-length, 
width Ly strips of the square lattice is quite different, even in the limit Ly ^ oo, from the 
behavior of the model on the square lattice. In the latter case, the thermodynamic limit is 
Lx oo, Ly oo, with limi^^oo Ly/ equal to a finite nonzero constant. For the strips, 
for any Ly no matter how large, the ferromagnet is critical only at T = 0, and as T ^ and 
^ — > oo, the strip acts as a one-dimensional system, since lirriL^^coLy/ L^ = 0. In contrast, 
for the Potts model on the square lattice, the phase transition occurs at finite temperature, 
at the known value i^c = ln(l + ^/q ). These studies of the thermodynamic behavior of the 
Potts model for general q on Ly x oo strips thus complement studies such as those on the 
approach to the thermodynamic limit of the Ising model on L^ x Ly rectangular regions, in 
which Lx and Ly both get large with a fixed finite ratio Ly/L^ [^, and finite-size scaling 
analyses |6^. These differences are also evident in the behavior of Bu] we have inferred above 



that as Ly oo, there are an infinite number of curves on Bu that cross each other at the 
ferromagnetic zero-temperature critical point, u = 0, so that the Fisher zeros become dense 
in the neighborhood of this point. This is quite different from the accumulation set of the 
Fisher zeros for the square lattice; although this is known exactly only for the Ising case, the 
existence of low-temperature expansions with a finite radius of convergence for the g-state 
Potts model is equivalent to the statement that the singular locus Bu does not pass through 
u = 0. 



In the case of the antiferromagnet, as we have shown |^T|, for q values that are only 
moderately above the value of g = 3 where the Potts antiferromagnet is critical on the square 
lattice, the ground state entropy of infinite-length, finite-width strips rapidly approaches its 
value for the square lattice. For the {L^ oo limit of the) Ly = 2 strip, this is given by 
5*0 = (l/2)fcBln(g^ — 3g + 3), which is nonzero for q > 2. The analytic expressions for the 
Ly = 3,4 cases are given in ^T|]. This can be understood because the ground state entropy 
is a disorder quantity and, for g > 3 is not associated with any large correlation length. 

2. < q < 2: Phase Transition for Antiferromagnet 

For the range < g < 2, our result for ac,+ in eq. ( |6.26| ) shows that B crosses the positive 
real a axis in the interval < a < 1, so that the Potts/random cluster antiferromagnet has 
a finite-temperature phase transition, at the temperature 

Tl,p = f , < g < 2 (6.45) 

fee In + V9^} 
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(where both J and the log are negative, yielding a positive Ti^p). For g = 1, it is understood 
that one takes n ^ oo first and then g — 1, i.e., that one uses the free energy As 
q decreases from 2 to 0, the phase transition temperature T^^p increases from to infinity. 
In the high- and low-temperature phases, the free energy is given by eq. ( |6.28|) and by 



/ = (l/2)lnAL,3, respectively. These results may be compared with the temperature Tp 
in eq. ( [4.14|) for the circuit graph. The same comments that we made in that case apply 
here; this result does not contradict the usual theorem that ID (and quasi-lD) spin systems 
with short-range interactions do not have any finite-temperature phase transition because 
the phase transition here is intrinsically connected with unphysical behavior of the model 
in the low-temperature phase, including negative specific heat, negative partition function, 
and non-existence of an n ^ oo limit for thermodynamic functions that is independent of 
boundary conditions. Indeed, the last pathology is obvious from the fact that for the n ^ oo 
limit of the ladder graph with open longitudinal boundary conditions, the free energy is given 
by eq. (|5.18|) for all temperatures, the singular locus B does not cross the positive a axis. 



and there is no such phase transition at finite temperature. 

Evidently, the temperature value at which the phase transition takes place in the 
Potts/random cluster antiferromagnet on the infinite- length limits of both the circuit graph 
and the cyclic and Mobius Ly = 2 strip graphs is determined by the respective formulas 
relating qc to a, eqs. ( [4.101 ) and (|6.25| ). From the point of view of Bq in the q plane, as 



we have discussed, we find, as a general feature, that in the antiferromagnetic case, as one 
increases T from to infinity, the value of qc{{G}) for a given family {G} decreases from its 
T = value to the origin, g = 0. Correspondingly, for the n ^ oo limit of a given family 
{G} with periodic (or twisted periodic) longitudinal boundary conditions, the antiferromag- 
net will exhibit a finite-temperature phase transition at a temperature 7{g},p for the range 
< g < g,({G}): 

3 T|G},p > for < g < g,({G}) . (6.46) 

Thus, for example, for the Ly = 3 square strip of the square lattice with cyclic or Mobius 
boundary conditions, for which we determined Bq for the T = antiferromagnet |^ and, in 



particular, gc(sg, Ly = 3, eye.) ~ 2.33654, it follows that the random cluster antiferromagnet 
has a finite-temperature phase transition for < g < qc{sq,Ly = 3, eye.). Just as we have 
discussed above, at special integer values g^ in the range < g < qc{{G}), it is understood 
that one takes the limit n — > oo first, and then g ^ g^ in calculating / = /g„ and Bu = {Bu)qn- 
Similarly, on the Ly = 2 cyclic and Mobius triangular lattice strips, where we found that 
qc{tri, Ly,eye.) = 3 for = 2 |^ and Ly = 3 [Q, it follows that the Potts/random 



cluster antiferromagnet has a finite-temperature phase transition for < g < 3. In all cases. 
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however, this transition involves unphysical aspects, among which is the non-existence of a 
unique n — > oo hmit that is independent of boundary conditions. 

I. Bu{{L}) for g > 4 

We next proceed to the shces of B in the plane defined by the temperature Boltzmann 
variable u, for given values of q, starting with large q. In the limit q oo, the locus i3„ is 
reduced to 0. This follows because for large q, there is only a single dominant Xj, namely 

Al,5 ~ + ^qv + 0{1) as q ^ oo . (6-47) 

Note that in this case, one gets the same result whether one takes g — > oo first and then 
n = 2m ^ oo, or n — > oo and then g — > oo, so that these limits commute as regards the 
determination of Bu- 
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FIG. 17. Zeros of Z{Lm, q, a) in the u = 1/a plane for q = 10 and m = 18 (n = 36). 
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FIG. 18. Zeros of Z{Lm,, q, a) in the u — 1/a plane for q = 3 and m — 24 {n — 48). 



In the large-g region we find that the locus Bu consists of two complex-conjugate curves 
that pass through m = at the angles (|6.23|) , hence intersecting at right angles and forming 
a distorted figure-8, together with a separate self-conjugate arc. Thus, Bu is comprised of 
two disconnected parts. For real q, the degeneracies in magnitude among leading terms 
are lA^^^a^^l = jA^^^s^ul at -u = and |Ai,_5^„| = |Ai^^6,u| at the point on the negative real axis 
where Tsu = 0. As a typical illustration of B^ for the large-g region, we show the complex- 



temperature zeros for g = 10 calculated for m = 18, i.e., n = 2m = 36, in Fig. |T^. It is 
instructive to compare this plot with the analogous plot for the open strip with g = 10 given 
above in Fig. |[ The self-conjugate arc is the same in the two plots, crossing the real u axis 
at u ~ —0.3954, where Tsu = and having endpoints at two of the zeros of the expression 
in the square root in Xlj, j = 5,6, which are identical to Xsj, j = 1,2 in eq. (|5.8|) . A 
notable difference between the locus B^ for the cychc or Mobius ladder and the analogous 
locus for the open square strip is that while the latter does not separate the u plane into 
different regions, the former does. Specifically, there are three regions: the physical PM 
phase occupying the interval < m < oo and its maximal analytic continuation, together 
with an Oi phase in the interior of the upper curve and its complex-conjugate phase 01- As 
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is evident in the figure, the density of zeros on the curves decreases strongly as u approaches 
the multiple point at m = 0. In the PM phase Xl,^,u is dominant, and so the reduced free 
energy is given by 

/ = ^lnAL,5 for uePM. (6.48) 

In the other phases only the magnitude \e-^\ can be determined unambiguously, and, with 
appropriate choices of branch cuts, we have 

|e^| = IAl.sT^^ for M G Oi, . (6.49) 

As q increases, Bu contracts toward u = 0, just as was true for the open strip. As q decreases, 
the point at which the self-conjugate arc crosses the negative real u axis (i.e. where Tsi2 = 0) 
moves toward the left, and the elongate toward the left. 



J. Bu for 3 



We next discuss the complex-temperature phase diagram for the case g = 3. An impor- 
tant conclusion that we shall draw from our studies of Bu for g = 3 and g = 4 (as well as 
the q = 2 case), building on our earlier comparative studies of complex-temperature phase 
diagrams for the ID and 2D Ising model with both spin 1/2 and higher spin s is that 

although ID and quasi-lD systems with short-range spin-spin interactions (infinite-length 
circuit or cyclic or Mobius strips) have qualitatively different physical thermodynamic prop- 
erties than the same systems in higher dimensions, the complex-temperature phase diagrams 
of these ID and quasi-lD systems can give insight into the corresponding phase diagrams of 
the model on lattices of dimensionality d = 2. Since no exact solution has been obtained for 
the Potts model in c? > 2 (except for the d = 2, q = 2 case), whereas we have exact solutions 
on infinite-length, finite-width strips, this means that one can use these as a tool to suggest 
properties of the complex-temperature properties of the Potts model in 2D (and perhaps in 
d > 2). 



The complex-temperature zeros of Z in the variable u are shown for g = 3 in Fig. |T8|. In 
addition to the (CTE)PM phase, which includes the intervals m > and u < —2 on the real 
u axis and the intervals 0.52 < \Im{u)\ < 1.7 and |/m(n)| > 2.2 on the imaginary axis, and 
extends outward to the circle at infinity, one also has several O phases. Among these are an 
Oi phase that contains the interval < Im{u) < 0.52 and an O2 phase which includes the 
interval 1.7 < Im{u) < 2.2, together with the complex conjugates of these phases, which are 
denoted 01 and The dominant terms in these phases are: Xl,5 in PM; Xl,3 in Oi, OJ; 
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and Xl,2 in O2, O^- Two phases that are self-conjugate and include intervals of the real axis 
are the O3 phase, containing the real interval — 1/2 < m < 0, in which A^^s is dominant; and 
the O4 phase, containing the real interval —2 <u< —1/2, in which A^^i is dominant. The 



point u = —2 here is the same as the point Uc{q) in eq. (|4.12| ) for the infinite-length limit 
of the circuit graph. There are also phases that have no support on the real axis. The locus 
Bu has several multiple points (in the technical terminology of algebraic geometry, meaning 
points where several branches of an algebraic curve intersect). Anticipating our results for 
other values of q, we find that the point on the negative real u axis where the PM phase, 
on the left, is contiguous with the O4 phase, on the right, is given by the same Uc as for the 
circuit graph, i.e., 

upM-oA{L}) = u,{{C}) = ^ . (6.50) 

q Z 

In a similar manner, we label the point on the negative real u axis where the O4 phase, on 
the left, is contiguous with the O3 phase, on the right, as uo^-o-i] as indicated, this has the 
value —1/2 for g = 3. The degeneracies in magnitude between leading terms \L,j,u at these 
multiple points are as follows (with appropriate conventions for branch cuts in square roots): 

|AL,3,n| = |Al,5,«| at M = (6.51) 

|Al,i,„| = |Al,5,^,| at u = uo^-os (6.52) 

|Al,i,u| = |Al,2,u| = |Al,5,u| at u = upM-Oi (6.53) 

|Al,i,J = |AL,3,n| = |Al,5,J at m ~ -0.44 ± 0.222 (6.54) 

\\L,j,u\ all equal at u = e^'^^l'^ . (6.55) 

Note that these points Ue are the same as for the open square strip discussed above. 

Motivated by our previous work ||5^j6^, we next explore relations between the exactly 



determined complex-temperature phase diagram for the Potts model on these strip graphs 
and on the square lattice. For g = 3 and higher integers, where the 2D g-state Potts model is 
not exactly solved, the complex-temperature phase diagram and associated Fisher zeros for 
various lattices have been studied in a number of works (e.g. ||66|- |7^, [^). One exact result 



concerning the complex-temperature phase boundary Ba for the g = 3 Potts model on the 
square lattice is that, as a result of the duality property ( ^.12| ) and the fact that the square 
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lattice is self-dual, Ba maps to itself when one replaces a with ad given by (|8.14|) ; hence 



the fact that the q = 3 Potts antiferromagnet on the square lattice has a zero-temperature 
critical point fl^, so that a = is on Ba, means that the dual of this point, namely a = —2, 



is also on Ba, or equivalently, u = —1/2 is on Bu Our exact results for the ID case with 
periodic boundary conditions, eq. ([4.11|) and for the infinite-length cyclic or Mobius ladder 
have the same feature, viz., that Bu contains the point u = —1/2: 

Bu3u = -1/2 for {C}, {L} and A.^ with g = 3 (6.56) 

where Agq denotes the square lattice. 

This interesting similarity of a feature of the complex-temperature phase boundary Bu 
leads us to investigate whether there are also other such connections. Our results suggest 
that the point u = —2 where for g = 3 the general formula ( |6.50|) shows that Bu{{G}) crosses 



the negative real axis for {G} = {C} and {L}, and the points u = e^^*'^/^ in eq. ( |6.55| ) where 
degeneracies in |Aj|'s and associated multiple points on B{{L}) occur, have analogues for 
Bu in the q = 3 Potts model on the square lattice. These conjectures could, in principle, 
be tested by calculations of Fisher zeros for g = 3 on finite square lattices, and these have 



been done; however, the considerable scatter of the zeros in the Re{a) < region |[66|- |69| , |7l 
renders it difficult to test the conjectures at present. One could also calculate the Potts 
model free energy and the boundary Bu on infinite-length strips of greater width, Ly > 3 
and check to see if for q = 3, the points u = —2 and u = e^'^'^'^l'^ are on the resultant locus 
Bu- For the analogous multiple points at m = ±z on Bu for the q = 2 Ising model we have 
done this (see below) and have found that these points are on this locus not only for Ly = 1 
and 2 but also for Ly = 3, just as they are for the full 2D square lattice. 

Let us comment further on the correspondences of special complex-temperature points 
for the present strip and for the square lattice. There are close relations with the proposed 
formulas (a— 1)^ = g for g > 2 and (a+1)^ = 4— g for 2 < g < 3 for the Potts ferromagnet 
and antiferromagnet, respectively. The second formula, (a + 1)^ = 4 — g (generalized to 
consider complex as well as physical temperatures) has the solutions a = — 1 ± ^y4 — q 
(which are duals of each other under the map a — = 1 + g/(a — 1). For g = 3, these are 
and —2. For g = 4, the two solutions coincide at a = u = —1. Both of these complex- 
temperature points, u = a"^ = — 2 for g = 3 and m = — 1 for g = 4, agree with eq. ( |6.50|) for 
the infinite cyclic or Mobius square strip. 

One of the earliest conjectures for the complex-temperature phase boundary Ba of the 
g = 3 Potts model on the square lattice was that it consisted of the square-lattice Potts 
model consists of the two circles |a — 1| = ^/g and |a + 1| = -^4 — g [^. From later 



studies of Fisher zeros, it was concluded that the Ba was not this simple ||67 H69H7l| . This 
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was further established combining the calculations of Fisher zeros with the analysis of low- 



temperature series expansions [[7T] , |f3[ , which showed the existence of complex-temperature 
singularities not on this locus, which were associated with prongs or cusps formed by the 
zeros. Nevertheless, if, indeed, the points a, a* = e^'^'^'^l'^ are on Ba-, (equivalently, Bu since 
the set of these points is the same under inversion) for the g = 3 square-lattice Potts model, 
as might be inferred from our exact results on the strips considered here, then this makes 
a very interesting connection with the old conjecture of Ref. |^6[, since the two circles are 
|a — 1| = and |a + l| = lforg = 3, and these intersect precisely in the two points 
a, a* = e^^'^*/^. We recall that if one uses self-dual boundary conditions, then one finds that 
the Fisher zeros lie nicely on the circle |^| = 1, where C, = [a — 1)/ ^/q, at least for Re{() > 



|69| and, moreover, in the g — oo limit, the complex-temperature phase boundary is \(\ = 1, 



whereC = (a-l)/yg 0. 



K. ^„ for 3 < g < 4 



As q increases in the real interval 3 < q < 4, the O4 phase contracts, as can be seen from 
the fact that the point upm-Oa inoves to the right, from —2 to —1, while the right-hand 
boundary at the point U04-O3 moves to the left, from —1/2 to —1; thus, at g = 4, these 
coincide: 

UPM-04 = U04-O3 = -1 for g = 4 . (6.57) 

In this interval 3 < g < 4, the degeneracies in magnitude of leading terms {XljI at upM-Oi 
in eq. ( p.53| ) and at ^04-03 in eq- ( |6.52|) continue to hold. For g = 4, all \\L,j\ are equal at 



u = —1. The confluence of upM-04 and M04-O3 at —1 and the equality of all lA^^jI's at this 
point reflect a special role for the value g = 4 for the complex-temperature phase diagram. 
(However, for the physical thermodynamics of these ID and quasi-lD systems, there is no 
qualitative change in the nature of the singularity at the zero-temperature critical point of 
the Potts ferromagnet at this value g = 4.) It may be recalled that the value g = 4 is also 
special, albeit in a different way, for the 2D Potts model in that for physical values of g below 
4 the Potts ferromagnet has a second-order phase transition while for g > 5 this transition 
is first order. 



L. ^„ for 2 < g < 3 

We next discuss the complex-temperature phase diagram as g decreases through real 
values from 3 to 2. From the point of view of this phase diagram, the limit g — 2 is singular. 
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since Bu is compact if g 7^ 2 but is noncompact for g = 2, passing through 1/u = 0. As q 
decreases through this range, the point upM-Oi moves leftward, approaching —00 as g — 2^". 
The point uo^-o-j, at which Bu crosses the negative real u axis separating the O4 phase on 
the left from the O3 phase on the right moves toward the right, from —1/2 to —0.453398.. 
(a root of the cubic equation + 2m + 1 = 0) as g decreases from 3 to 2. In Fig. [2^ we 
show a plot of complex-temperature zeros for g = 2.5. One can observe how the intersection 
points which occurred at m = e^^*'^/^ at g = 3 have shifted outward from the real axis and 
toward the right. When g decreases all the way to 2, these intersection points reach ±z (see 
below). For g = 2.5, the point upM-Oi = while the crossing at m = uq^-o-j, is clearly 
visible, near to its g = 2 limiting location at u ~ —0.4534. 




FIG. 19. Zeros of Z{L„n q, a) in the u = 1/a plane for q = 2.5 and m = 24 (i.e., n = 48). 
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FIG. 20. Same as Fig. |9|for g = 2. 



M. Bu for g = 2 



Here we again encounter noncommutativity in the definition of tfie free energy and the 
locus Bu- We first discuss fnq and {Bu)nq, obtained by setting q = 2 and then taking n ^ oo. 
For q = 2, besides the g-independent A^^i = (a — 1)^, we have 



Al,2 = (a - l)(a+ 1) 



(6.58) 



A 



L,(5,6) 



o(« + l) 



Al,3 = a{a- l)(a + 1) 



Al,4 = (a - 1) 



+ 1± (a^ -Aa-^ + lOa^ - 4a + 1 



1/2 n 



(6.59) 
(6.60) 

(6.61) 
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so that for this value of q, = \l,4:- Further, cl,i = —cl,^ = —1 so that the {Xl,i)^ and 
(•^^,4)™ terms cancel each other in Z, which reduces to 

Z(L^,g = 2,a)= {K.T ■ (6-62) 

j=2,3,5,6 

Since each of the A/s contributing to Z{Lm, q = 2,a) has an (a + 1) factor, it follows that 

Z{Lm, q = 2,a) = 2(a + I)*" x polyn. (6.63) 
The result (|6.63|) implies that Z{L, q = 2,a = —1, m) = 0, which is implied more generally by 



our earlier Theorem 6 of Ref. |5^. This theorem states that for lattices with odd coordination 
number, the (zero-field) partition function of the Ising model vanishes at a = — 1. This point 
may or may not be on {Bu)nq', for example, for the honeycomb lattice, it is, while for the 
Archimedean 3-12 lattice it is not [Q. In the present case we will find that it is not on 

In Fig. |2D| we show a plot of Fisher zeros for the q = 2 case. Because of the general 
fact that {Bu)nq passes through u = and the special inversion symmetry (|2.4|), ( p.l6| ) that 
holds for g = 2, it follows that {Bu)nq passes through a = also. The complex-temperature 
phase diagram consists of six phases, bounded by the locus Bu which is the continuous 
accumulation set of the Fisher zeros. The phase diagram in the u plane consists of (i) the 
complex-temperature extension of the Z2-symmetric, paramagnetic phase, PM, including the 
real axis m > 0, together with five O phases: (ii) Oi, including the interval on the imaginary 
u axis from u = to u = i, and its complex-conjugate, (iii) 01; (iv) O2, including the interval 
on the imaginary axis from u = i to u = zoo, and (v) its complex conjugate, O2', and (vi) O3, 
including the negative real axis, m < 0. The PM and O3 phases map to themselves under 
the inversion u —>■ 1/u = a, while Oi ^ O2 and O2 ^ 01 under this inversion. As in our 



previous work |S^JS^ we shall henceforth suppress the qualifier (CTE) and refer to the PM 
phase simply as the PM phase. In this phase, A^^s is the dominant term in ( |6.1| ) so that 
reduced free energy is given by / = (l/2)lnA2,^5 in the PM phase, as in ( |6.48| ). In other 
regions that are not analytically connected with the PM region, only |e-^| can be determined 
unambiguously: 

|e^| = |AL,5r/' for we 03,0* (6.64) 



Ai,3|^/2 for ueOi,Ol 



(6.65) 



and 



60 



e 



/| = \Xl,2\^/^ for ue 02,0* . (6.66) 



where here / = The locus Bu has multiple points (in the algebraic geometry sense) 

at M = and u = ±i and at complex infinity, i.e., ata = l/u = corresponding to the 
zero-temperature ferromagnetic and antiferromagnetic Ising critical points. We have also 
calculated Z{G,q,v) exactly for the Ising model on the cyclic Ly = 3 strip and have found 
that Bu again has multiple points at u = ±z (as well as at u = and 1/u = 0). A possible 
inference would be that this is true for cyclic or Mobius strips of the square graph for all 
finite values of Ly > 2. Given our exact results for Ly = 2, 3, we see once again a very 
interesting connection with the complex-temperature phase diagram of the same model - in 
this case, the Ising model - on the two-dimensional square lattice, for which Ba consists of 
the Fisher circles |m ± 1| = V^, which intersect precisely in the points u = ±i. 

If one takes n ^ oo first and then g — > 2, the resulting locus {Bu)qn differs from (Bu)nq in 



several respects. First, {Bu)qn does not satisfy the inversion symmetry ( p.l6|) . Second, while 
^L,5,u is dominant on the negative u axis in the vicinity of the origin u = 0, \l,i,u (equal, in 
the g — > 2 limit, to A^, 4^^) becomes dominant for u < —0.454 and similarly for radial paths 
emanating outward from the origin in the upper and lower Re{u) < half-plane. This gives 
rise to another region boundary. (Recall that the contributions of these A's cancelled if one 
took q = 2 first, so that they did not affect Z or the locus {B)nq in the n — >■ 00 limit.) The 
absence of the inversion symmetry in {Bu)qn is clear since a region boundary on this locus 
passes through u ~ —0.454 but not the inverse of this point. This is, then, an example of 
the noncommutativity {Bu)qn 7^ {Bu)nq for a case where both of these loci are nontrivial. 
Finally, one can also discuss Bu for negative real q and for complex q, but we shall forgo this. 



VII. SUMMARY AND CONCLUSIONS 



In summary, we have calculated exact closed-form expressions for the Potts 
model/random cluster partition function for general q and temperature T, or equivalently, 
the Whitney/Tutte polynomial for the open, cyclic, and Mobius square strips (ladder graphs) 
of width Ly = 2 and arbitrary length L^. Taking the limit L^ 00, we have determined 
the free energy / (and \e^ in unphysical phases) and the continuous locus B where the free 
energy is singular, which arises as the continuous accumulation set of the partition function 
zeros in the space of the variables q and u. The divergences in the correlation length and 
susceptibility of the Potts ferromagnet at its zero-temperature critical point were shown to 
be more rapidly approached as the strip width increases, and the physical reason for this was 
given. Our comparison of different strip widths suggests the inference that for infinite-length. 
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width Ly cyclic (or Mobius) strip graphs, as Lj, — > oo, an infinite number of curves on Bu 
pass through the point m = and the Fisher zeros become dense in the neighborhood of this 
point. It was shown that the Potts/random cluster antiferromagnet on both the infinite- 
length circuit graph and ladder graph with cyclic or Mobius boundary conditions exhibits a 
phase transition at finite temperature if < g < 2, but with unphysical properties. We dis- 
cussed a subtlety in the definition of the free energy of the random cluster model due to the 
noncommutativity at certain special values g^: lim„^oo liniq^g^ Z^/^ ^ limg^^^ lim^^oo 2'^/". 
Several generalizations of results for the T = limit of the Potts antiferromagnet (chromatic 
polynomials) were presented. Among these is the general form (p.l8|) for the partition func- 
tion of recursive strip graphs (and its further generalization to the case of nonzero external 
field, (|2.24|) ). The analysis of ||25| for the singular locus Bq in the case of chromatic polyno- 
mials was generalized to the present case of the full temperature-dependent Potts partition 
function. The dependence of the locus S as a function of the longitudinal boundary con- 
ditions was studied, and it was shown that this locus is the same for the strips considered 
here with cyclic and twisted cyclic (Mobius) longitudinal boundary conditions. For the Potts 
antiferromagnet, it was found that as the temperature increases from to infinity, the sin- 
gular locus Bq contracts in to the origin, g = 0. For the Potts ferromagnet, this locus was 
found not to cross the positive q axis, in contrast to the antiferromagnetic case, where, for 
the cyclic strip, it does. In both the antiferromagnet and ferromagnet cases, for the strips 
studied here, Bq passes (does not pass) through g = if one uses periodic (free) longitudinal 
boundary conditions. Generalizing our previous result for chromatic polynomials, we found 
that for the strips with periodic longitudinal boundary conditions, Bq encloses regions in the 
g plane for values of a where it is nontrivial (i.e., a 7^ 1). Several advantages of periodic, as 
opposed to free, longitudinal boundary conditions were noted, including the fact that with 
such periodic longitudinal boundary conditions, the locus Bu passes through m = 0, in 1-1 
correspondence with the zero-temperature critical point of the Potts ferromagnet. Finally, 
certain properties of the complex-temperature phase diagrams and loci Bu for these infinite- 
length, finite-width strips were shown to be the same as known properties of the model on 
the square lattice, including the multiple points at m = ±2 for g = 2, which also occur in 
the exactly solved square-lattice Ising model and the point u = —1/2 for g = 3, which is 
known to lie on B^ for the square lattice since it is dual to the zero-temperature critical 
point at a = 0. This shows that exact solutions on infinite- length strips could provide a way 
of generating plausible conjectures for complex-temperature properties of the Potts model 
on two-dimensional lattices and some conjectures were made. 
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VIII. APPENDIX 

A. Connection Between Potts Model Partition Function and Tutte Polynomial 

The Potts model partition function Z{G, q, v) is related to the Tutte polynomial T(G, x, y) 
as follows. The graph G has vertex set V and edge set E, denoted G = {V, E). A spanning 
subgraph G' is defined as a subgraph that has the same vertex set and a subset of the edge 
set: G' = {V,E') with E' C E. The Tutte polynomial of G, T{G,x,y), is then given by 0- 

i 

nG,x,y)= J](a;-l)'=(G')-fc(G)(^_i)c(G') ^g^) 

G'CG 

where k{G'), e(G'), and n{G') = n{G) denote the number of components, edges, and vertices 
of G", and 

c{G') = e(G") + k{G') - n{G') (8.2) 

is the number of independent circuits in G' (sometimes called the co-rank of G'). Note that 
the first factor can also be written as (x — \y'y^)~^'y^') ^ where 

r{G) = n{G) - k{G) (8.3) 

is called the rank of G. The graphs G that we consider here are connected, so that k{G) = 1. 
Now let 

X = 1 + - (8.4) 

V 



^Accordingly, the U.S. government retains a non-exclusive royalty-free license to publish or reproduce the 
published form of this contribution or to allow others to do so for U.S. government purposes. 
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and 

y = a = v + l (8.5) 
so that g = (x — — 1) = (x — l)v. Then 

Z{G, q, v) = {x- {y - \)<^^T{G, x, y) . (8.6) 
There is also a connection with the Whitney rank polynomial, R{G,^,ri), defined as P,p!^ 

R{G,^,v)=J2c'^'k'^'''^ (8.7) 

G'CG 

where the sum is again over spanning subgraphs G' of G. Then 

T(G, X, y) = {x- l)'-(^)/?(G, e = (x - l)-\ = 2/ - 1) (8.8) 

and 

Z{G, q, v) = g"(^)i?(G, ^ = -,r^ = v). (8.9) 

Note that the chromatic polynomial is a special case of the Tutte polynomial: 

P(G, q) = q^i(^)(-lfiG)+n{G)rj.^Q^ 2; ^ 1 _ ^ ^ 0) (8.10) 



(recall eq. ([LTOD ). 



From the representation ( ^.61) and the duality property of the Tutte polynomial 

T{G,x,y)=T{G\y,x) (8.11) 
where G* is the dual graph corresponding to G, it follows that 

Z(G, g, t;) = t;^(^)g-=(«)Z(G*, g, t;,) (8.12) 
where G* denotes the graph that is dual to G and Vd is the dual image of v: 

Vd = - (8.13) 

V 

or equivalently, in terms of the variable a, 

a — 1 + g , , 

ad = . 8.14 

a — 1 

Corresponding to the form ( |2.18| ) we find that the Tutte polynomial for recursively defined 
graphs comprised of m repetitions of some subgraph has the form 

T{Gm,x,y) = Y,^t,g,Mt,gT (8-15) 
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B. Square Strip with Free Longitudinal Boundary Conditions 



The generating function representation for the Tutte polynomial for the open square strip 



We have 



where 



m=0 



TT{S,x,y]z) 



AfriS^x^y; z) 
VTiS,x,y;z) 



(8.16) 



(8.17) 



JVt{S, X, y; z) = AT,sfi + At,s,iz = {y + x + x + x ) - yx z 



and 



Vt^S, x,y, z) = 1 — {y + 1 + X + x^)z + yx'^z^ 



with 



3=1 



Ar,s,(i,2) = \ {l + y + x + x^)±(y'^ + 2y{l + x - x^) + {x^ + x + I) 



l/2n 



(8.19) 



^.20) 



The corresponding closed-form expression is given by the general formula from [^], as applied 
to Tutte, rather than chromatic, polynomials, namely 



T(Sm,x,y) 



At,s,o^t,s,i + At, 



5,1 



^T,S,1 ~ ^T,S,2 



A-T,S,0^T,S,2 + At,S,1 



(At,5,2)™ . (8.21) 



An alternative expression for T that explicitly shows that it is a symmetric function of 
the Asj, j = 1,2, is 



T{Sm, x,y) = ^{y + x + x^ + x') (At.s.i)'" + (At,5,2)' 



+ 



y^ + y + 2yx + 2yx^ + x + 2x^ + Sx"* + 2x'' - yx^ + x^ 



i^T,S,l)"^ — (At,5,2) 



At.5,1 ~ A 



T,5,2 



(8.22) 
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C. Cyclic and Mobius Square Strips 



We write the Tutte polynomials for the cyclic and Mobius square strips Lm and MLm as 



and 



T{MLm,X,y) = ^CT,ML,ji>^T,ML,jy 



where it is convenient to extract a common factor from the coefficients: 



""^'^'^ G = L,ML . 



X — 1 



(8.23) 



(8.24) 



(8.25) 



Of course, although the individual terms contributing to the Tutte polynomial are thus 
rational functions of x rather than polynomials in x, the full Tutte polynomial is a polynomial 
in both X and y. We have 



^T,ML,j — ^T,L,j , j — 1, 6 



(8.26) 



At,l,i — 1 



(8.27) 



and 



^T,L,2 — X 



A 



T,L,(3,4) 



X + y + 2 ± (^{x - y f + 4:{x + y + 1) 



1/2 



(8.28) 



1.29) 



At,l,5 — Ar,s,i 5 At,l,6 — At,5,2 • (8.30) 
Our result for T{G, x,y), G = L, ML agrees with a recursion relation given in p2[ (see also 



cr,L,i = [{x -l)iy- I)? - 3{x - l)(y - l) + l 



(8.311 



ct,l,2 = cr,L,3 = ct,l,4 = xy - X - y 



(8.32) 
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Cr,L,5 — Ct,L,6 — 1 

ct,ml,i = ~1 
ct,ml,2 = -xy + x + y 
ct,ml,3 = cmla = xy - X - y 



8.33) 
8.34) 
8.35) 
8.36) 
8.37) 



Ct,ML,5 — CT,MLfi — 1 ■ 

We note that AT,L,3Ar,L,4 = xy and Xt,l,5>^t,l,6 = x'^y. 

D. Special Values of Tutte Polynomials for Square Strips 

For a given graph G = (V, E), at certain special values of the arguments x and y, the Tutte 
polynomial T{G,x,y) yields quantities of basic graph-theoretic interest ll^, citewu77- 



||79| . We recall some definitions: a spanning subgraph G' = (V, E') of G is a graph with the 
same vertex set V and a subset of the edge set, E' C E. Furthermore, a tree is a graph 
with no cycles, and a forest is a graph containing one or more trees. Then the number of 
spanning trees of G, Nst{G), is 

Nst{G) = T{G, 1,1) , (8.38) 
the number of spanning forests of G, Nsf{G), is 

Nsf{G) = T{G,2,1) , (8.39) 
the number of connected spanning subgraphs of G, Ncssg{G), is 

Ncssg{G) =T{G, 1,2) , (8.40) 
and the number of spanning subgraphs of G, Nssg{G), is 

Nssg{G) = T{G,2,2) . (8.41) 

Clearly, NssGiG) — Ncssg{G) is the number of disconnected spanning subgraphs of G and 
Ncssg{G) — Nst{G) is the number of connected spanning subgraphs of G that contain one 
or more cycles. One thus has the inequality 
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NssciG) > Ncssg{G) > Nst{G) i.e., T{G, 2, 2) > T{G, 1, 2) > T{G, 1, 1) 



(8.42) 



Also, clearly 



NssaiG) > NsFiG) i.e., T{G, 2,2) > T{G, 2,1) 



(8.43) 



and 



NsFiG) > Nst{G) I.e., T(G, 2, 1) > r(G, 1, 1) . 



(8.44) 



The set of spanning forests differs from the set of connected spanning subgraphs by the 
removal of the condition that the subgraph is connected but the imposition of the condition 
that the subgraph have no cycles, and hence there is no general inequality between Nsf{G) 
and Ncssg{G). 

We recall the results for tree graphs T„ and circuit graphs C„: T(Tn, 1, 1) = T(Tn, 1, 2) = 
1, T(T„, 2, 1) = T(T„, 2, 2) = 2--\ T(C„, 1, 1) = n, T{Cn, 2, 1) = 2" - 1, T(C„, 1, 2) = n + 1, 
and T{Cn, 2, 2) = 2". 

For the open square strip we find 



Nst{S„ 



(2 + V3 )" + (2 - V3 



+ 



2VS 



(2 + 73 )"^ - (2 - V3 



^.45) 



15 



NsFiSm) = - (2(2 + VS))'- + (2(2 - v^))' 



+ 



13 

7! 



(2(2 + 73))"^ - (2(2 - v^))" 



(8.46) 



n 



and 



5 + V17 \™ /5 - Vl7 ^ ™- 

+ 



21 
2717 



5 + V17 



^.47) 



Nssg{S^) = 2 



3m+4 



(8.48) 



That eqs. (|8.45| )- (|8.471) yield integers follows from the theorem on symmetric polynomial 
functions of roots of an algebraic equation, as discussed in [^|. 
With the definition 



Vg 



-1 ifG = L 
-1 if G = ML 



(8.49) 
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our calculations of the Tutte polynomials for the cyclic strip and the Mobius strip MLr, 
yield 



NsriGm) = mj-r/c; + ^ [(2 + Vs)" + (2 - Vs)™] } , = L 



-'mi -^Lfyi 



NsFiGm) = Vg{1 - 2' 



r/S + Vl7 \™ _^ /5 - Vl7 ^ ™' 



+ (^2(2 + 73))"+ (2{2 



} Gm — MLm 



(8.51) 



NcssciLm) = Ncssg{MLJ - 2m - 1 = -(m + 2) + 



(m + 2) 



m 



2V17 



17 \ ' 



and 



(8.52) 



Nssg{L^) = NssciMLr 



■)3m 



(8.53) 



The results for the spanning trees for the cyclic and Mobius strips are known (for higher- 
dimensions, see |T9|); we are not aware of the other quantities having been published. 

Several comments are in order. Since T{Gm,x,y) grows exponentially as m — > oo for 
the families Gm = Sm, L^, and MLm for {x,y) = (1,1), (2,1), (1,2), and (2,2) (as well as 
for Gm = Gm and Tm for {x,y) = (2,1) and (2,2)), it is natural to define corresponding 
constants 



Zset{{G})= lim n{G)-HnN,{G) 

n{g)^oo 



set = ST, SF, CSSG, SSG 



(8.54) 



where, as above, the symbol {G} denotes the limit of the graph family G as n{G) oo (and 
the z here should not be confused with the auxiliary expansion variable in the generating 
function ( ^.16| ) or the Potts partition function Z{G,q,v).) The general inequalities ( |8.42| ), 
( ^.43| ), and ( |8.44| ) imply that, for a given {G}, 



ZSSG > ZcsSG > ZsT 



1.55) 



ZSSG > ZSF 



^.56) 



and 
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ZSF > ZST 



(8.57) 



We find that for both the hne {Ly = 1) and the Ly = 2 square strip, the quantity 
Zset{{G}) is independent of whether the longitudinal boundary conditions are free, periodic, 
or Mobius: 

zst{{G}) = - ln(2 + v^) ^ 0.658479 for G = S,L, ML (8.58) 
2 

zsf{{G}) = ^\n[2{2 + V3)]- 1.00505 for G = S,L,ML (8.59) 



X / 5 ~l~ /IT 

zcssg{{G}) = - In ( \ ) ~ 0.758832 for G = S,L, ML (8.60) 



and 



2 V 2 



255g({G'}) = ^ In 2 ~ 1.03972 for G = S,L,ML (8.61) 



(where the result for Zst can be extracted from ^^^). 

E. Tutte Polynomials for Dual Graphs 



Since the Tutte polynomial satisfies the duality relation (^.11), our calculations of the 



Tutte polynomials T{Gm, x, y) for the open, cyclic, and Mobius square strips, G^ = Sm, L^, 
and MLm, also yield the corresponding results for the duals of these graphs. The dual of 
the square Ly = 2 strip with m + 1 squares, i.e., length L^ = m + 1 edges, (5'^)*, can be 
described as follows: for m > 1, consider a line of m + 1 vertices, with successive vertices f « 
and fj+i connected to each other by an edge e^; the vertices on this line are connected to a 
single external vertex by double edges, except for the first and last vertices on the line, each 
of which is connected to the external vertex via three edges. This is (Sm)* for m > 1. Note 
that this is a multigraph, since a (proper) graph is normally defined not to have loops or 
multiple edges. For the case m = 0, i.e., a single square, the dual is (^i)* = TL4, where in 
the mathematical literature, TLi, denoted "thick link", is the multigraph consisting of two 
vertices connected by i edges. Our results for T{Sm,x,y) in eq. ( |8.21| ) thus give the Tutte 
polynomial for this dual graph as T{{Sm)*,x,y) = T{Sm,y,x). 

For the dual of the cyclic strip graph, (Lm)*, we recall a definition from graph theory: 
given two graphs G and H, the "join" G + H is the graph obtained by connected each vertex 
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of G to each vertex of H with edges. We also recall the notation Kp for the complement of 
Kp, i.e. the graph consisting of p vertices with no edges. Then for m > 3, 



{LmT — K2 + C„ 



(8.62) 



where is the circuit graph with m vertices. Hence our results for T{Lm, x, y) in eqs. ( 6.1 
with ( |6.4|) - (|6?l0| ) also determine T{{Lm)* ,x,y) = T{Lm,y,x). 
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